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Abstract: In this brief review we will discuss how a well motivated particle the- 
^ ory beyond the eletroweak Standard Model provides ingredients and conditions for a 

O ! successful inflation. We will mainly focus on a low energy supersymmetric Standard 

vn . 

en ■ Model which provides plenty of scalars. In particular, these scalars span a multi- 

rn 

O \ Model charges. The inflationary predictions which matches the current observations 

O ■ are robust due to the fact that inflation occurs within our own gauge sector where 

the couplings are well known. We further argue that based on our current under- 
standings if there exists a string landscape of multiple vacua, then it is very natural 
that the last phase of inflation would be driven by one of the many supersymmetric 
Standard Model modulii. Only such a graceful exit from inflation would provide 
hot thermal Standard Model baryons, cold dark matter, conditions for baryogenesis 
and foremost the seed density perturbations for the cosmic microwave background 
radiation in just one package. Furthermore we will also discuss how some of the 
ingredients of inflation can be tested already by the LHC 
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1. Introduction 

Primordial inflation as a paradigm has met with glorious successes over the past 26 
years since its advent Jl], ^ ^, |^ (for a review see 0). The virtues of inflation 
lies due to the fact that an accelerated expansion of the universe can give rise to 
a homogeneous and an isotropic universe on very large scales with a flat spatial 
geometry. There are very important observational consequences, such as nearly scale 
invariant tiny perturbations imparted to the cosmic microwave background radiation 
and its spectrum [|] (for a review see 0) has met with an unprecedented success 
with the observations based on Cosmic Microwave Background (CMB) radiation, see 



the recent data from WMAP |T0|, In 



It has been known for a while that inflation can be driven by a dynamical scalar 
field known as the inflaton, an order parameter, which could either be fundamental or 
composite. Particularly if the inflaton rolls very slowly on a sufficiently flat potential 
(such that the potential energy density dominates over the kinetic term) then it is 
possible to show: (1) a considerably large number of e-foldings, (2) during inflation 
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metric perturbations are generated with an amplitude: Vr ~ 1.91 x 10 ^A;" ^, with 
a spectral tilt lying within: 0.92 < n < 1.0 within 2a (CMB+SDSS) (3) running of 
the spectral tih is negligible -0.012 < dus/dlnk < 0.001 within 2a (CMB+SDSS 
for tensor to scalar ratio negligible) and, (4) the tensor to scalar ratio, r < 0.55 at 
95% c.l. details can be found in Ref. [[10|] . 

However, inspite of all these achievements, there are some fundamental issues 
where we are lacking proper understanding, such as: 

• What is the origin of inflation ? 

• What is the inflaton ? 

• What are the fundamental interactions of an inflaton ? 

• Where does the inflaton energy goes ? 

• How can we test the inflaton in a laboratory ? 

Inspite of many attempts there has been no single good candidate for an inflaton 
which comes naturally out of a well motivated theory of particle physics (for a review 
on models of inflation, see |]12|). One always relies on scalar fields which are absolute 
gauge singlets possibly residing in some hidden sector or secluded sector with a small 
coupling to the SM gauge group. By definition an absolute gauge singlet does not 
carry any charge what so-ever be the case. Therefore the masses, couplings and 
interactions are not generally tied to any fundamental theory or any symmetry. Such 
gauge singlets are used ubiquitously by model builders to obtain a desired potential 
and interactions at a free will in order to explain the current CMB data. In this 
respect the inflaton is just a phenomenological answer to inflation and its observed 
consequences. What has been ignored so far is the fact that if inflation explains 
the CMB anisotropy then it must also reheat the universe with the SM degrees of 
freedom ^. The decay of inflaton should also provide thermal conditions |]18[ to 



^An absolute gauge singlet inflaton can couple to all the species and sectors alike. It does not 
discriminate between the SM degrees of freedom to that of any other degrees of freedom. By 
definition an absolute gauge singlets can couple to any other gauge singlets, there is no symmetry 
which prohibits doing so. The most popular slow roll inflationary models such as chaotic inflation 
hybrid inflation , new inflation [||, || , assisted inflation 1 14 , ^ , without slow roll inflation jl6j , 



K-inflation ||T^, etc. all rely on absolute gauge singlets, see for a review [|2|. In this respect they all 
fail to explain why the couplings, masses and interactions are chosen so to explain the CMB data. 
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create matter- anti- matter asymmetry and thermal/non-thermal generation of cold 
dark matter, for a review see [1^, ^. 

On the other hand string theory comes up with a plethora of absolute gauge 
singlet modulii, which mainly arise in the gravitational sector upon various compact- 
ifications . There are numerous attempts to embed inflation within string theory, 

. However they all tend to concentrate on ex- 



for a review see |2^, see also p5| . 



plaining the CMB data, but with very little to do with an aftermath of inflation 
such as how the infiaton couples to the SM degrees of freedom ?, how the dark matter 
is created ?, what are the relativistic degrees of freedom after the end of inflation 
?, etc. The uncertainties prevail mainly due to two facts: (a) It is not clear how to 
construct SM like gauge group. There are various constructions so far, but all suffer 
through various problems such as quantum instability, moduli de-stabilization, extra 
U(l)'s, etc. [^. (b) Inflation again generically driven by absolute gauge singlets. 



such as the closed string modulii in the case of a race-track inflation or open 



string moduli [29, 30], for a recent discussion on open string moduli inflation in a 



warped throat, see In the case of an open string moduli the infiaton need not be 
completely an absolute gauge singlet, but so far a realistic model of brane-anti-brane 



inflation lacks the SM embedding |32 



In all these cases the infiaton candidate can couple to closed string degrees of 
freedom and the open string degrees of freedom alike. Although there will be some 
preference to decay channels for an open string moduli as an infiaton, it can couple to 
an appropriate gauge theory with appropriate gauge couplings |^ . However there is 
no obvious construction which connects open string moduli as an infiaton to that of 
the SM sector. Moreover because of a hierarchy between the four dimensional Planck 
scale, string scale and the compactification scale, the 6 dimensional compactified 
volume is so large that a kinematical phase factor argument allows the infiaton 
energy density to get lost in the bulk of the volume which contains mainly the 



closed string degrees of freedom ||30| , |33| . Therefore, there is no preference why 
the infiaton would decay primarily into the SM baryons. One would have to make 
Moreover one cannot predict definitely certain issues such as radiative corrections and stability of 
the potential, because the couplings and masses are all put by hand. 

^The Big Bang Nucleosynthesis (BBN) demands that when the universe was nearly 1 MeV, most 
of the energy density of the universe must lie in the SM relativistic degrees of freedom . It is 
therefore pertinent that the universe must have created a thermal bath of strictly speaking SM 
degrees of freedom before the BBN. If there were other degrees of freedom which were to couple 
very weakly to the SM then it would provide severe experimental constraints |22||. 
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additional assumptions regarding various sectors the inflaton couples to or where 
exactly inflation occurs ^. The string motivated models have thus failed to explain 
the aftermath of inflation, particularly why the universe is filled with the SM baryons 
at the time of BBN. More work is required to understand this issue, until then 
one could at best treat stringy inflation as a nice playground to test various ideas 
concerning CMB data alone. 

There are also some attempts to understand the CMB data without invoking 
inflation such as in the case of a bouncing cosmology ^7\, . Some of the key 
ideas are stringy in origin (based on t-duality, see ||3^), especially in the context of 
string gas cosmology of Brandenberger and Vafa However, there are various 

caveats towards a full understanding of a thermodynamics of a gas of strings in a 



cosmological background pO[. There are other attempts to obtain inflation such as 



in the case of a gas of D-branes or string moduli in conjunction with string 



gas or in the case of Hagedorn strings during a bounce Moreover they all 
fail to address one crucial issue: how to obtain the SM degrees of freedom at the end 
of the day? 

An an important message from the above discussions is following: 

// any model of inflation or any alternatives of inflation provides the right CMB pre- 
dictions then there must be a way to excite the SM baryons with a precision which 
can match the BBN data without making further assumptions. 

In order to facilitate this we need a guiding principle, some kind of symmetry ar- 
gument, which will help us in constructing a successful model of inflation from a 
bottom-up approach. 

Very recently an interesting construction of inflation has been suggested within a 
minimal extension of the SM such as in the case of Minimal Supersymmetric Standard 
Model (MSSM) |g, ||, |7HS|, g (for a review on MSSM, see [0, ^ and 



its embedding in supergravity (SUGRA), see ||5J, IS^j, for cosmological aspects of 



MSSM, see |]T9|, Within MSSM there are many scalars, which span into a moduli 

^In a warped throat geometry, motivated by stabilizing the dilaton, complex structure moduli! 
and the kahler modulii, the inflaton sector and the SM sector are kept in distinctive throats. The 
SM throat typically addresses the hierarchy issue while the inflationary throat is situated relatively 
at higher energies. The two throats talk to each other through closed string sector. There are many 
issues which plague this geometry, such as excess Kaluza Klein graviton production |3|] , and 
transferring the energy density from closed string sector to the SM throat |35 . There are various 
subtleties as during inflation the SM modulii are displaced away from their minimum which brings 
in various cosmological uncertainties which might as well hamper baryogenesis and BBN ||35|| . 
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space of gauge invariant F-and D-flat directions, which carry SM charges usually 
baryon and/or lepton number. These modulii have an enhanced symmetry point 
near the origin (at a VEV defined by zero). Away from the origin these modulii 
break wholly or partly the SM gauge symmetry depending on the fiat direction. 
But such a spontaneous breaking of charge and color in the early universe is not 
considered to be dangerous, provided they all settle down to their minimum before 
the electroweak phase transition. 

In the first paper it was pointed out that the MSSM provides all the neces- 



sary key ingredients for a successful infiation. Out of nearly 300 fiat directions pS 



there are only 2 directions which can support infiation with a graceful exit. Other 
directions do not have a graceful exit |^, nor do they have a slow roll phase 
of infiation. This is a remarkable result which puts forward the two infiaton candi- 
dates as LLe (which carries the lepton number) and, udd (which carries the baryon 
number). The L corresponds to the slepton (SUSY partner of SM lepton) doublet, 
while e corresponds to the right handed selectron (SUSY partner of electron). The 
u, d, d correspond to the right handed squarks (SUSY partner of quarks). These 
two directions can support eternal and slow roll infiation with a right amplitude of 
temperature anisotropy Vr ~ 5// ~ 1.91 x 10~^, an observable tilt in the power 
spectrum with a range 0.92 < Ug < 1.0, where the lower limit is saturated for a sad- 
dle point infiation, the upper limit corresponds to having a slow roll infiation away 
from the saddle point [0, The running of the tilt in the spectrum is negligible 
and the tensor to scalar ratio is also observationally negligible. Moreover there is no 
production of cosmic strings, or non-MSSM degrees of freedom ^ , and there is 
no large non-Gaussianity |[57||. 



Furthermore the infiaton reheating is very well understood as the infiatons carry 
the SM charges, they naturally decay only to the SM degrees of freedom, besides 
creating the lightest supersymmetric particle (LSP) as a candidate for the cold dark 



matter With an R-parity LSP's are absolutely stable and can be a SUSY 

cold dark matter candidate 0. In Ref. we studied the parameter space of 
infiation and the neutralino type dark matter produced thermally. Remarkably we 
found an interesting overlap between infiationary parameters within mo, mi/2 (soft 
masses and gaugino masses respectively) plane in (m-SUGRA) setup. This provides 
another hint that if the desert is filled by MSSM from the electroweak scale to 
the grand unified scale, then the parameters of an MSSM infiaton will provide not 
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only a successful inflation which matches the current observational data, but the 
same parameter space is also responsible for generating the observed abundance of 
thermally created cold dark matter . 

Both the infiatons, LLe, udd, masses are tied with a low energy SUSY required 
to address the hierarchy between the Planck and the electroweak scale. The mass 
range ~ 0{TeV) is adequate to explain the CMB anisotropy. An arbitrary increase 
in the scalar masses will ruin the CMB predictions ||4^, ^ . Not only this, we would 



also be able to put constraints on the slepton and squark masses from the LHC [4S 



This will further constrain the MSSM inflation model to an unprecedented level. 
In this regard both LLe and udd are testable through CMB and the LHC. This is 
the most exciting expectation that for the first time the LHC will be able to rule 
out a model of inflation completely. For instance if the slepton and squark masses 
are beyond ^ 10 TeV, then this model of inflation is ruled out completely. For 
such a mass range neither LLe nor udd can generate the right amplitude of density 
perturbations ^ . 

Furthermore this is the only known model of inflation where the A/cobe (the last 
number of e-foldings required to explain the CMB data) is determined fully from the 
fundamental interactions. The A/cobe is primarily determined by thermal history of 
the universe. Fortunately enough, in MSSM based inflationary model, the reheating 
and thermalization temperature can be estimated rather accurately due to known 



SM couplings [||, H, |6|, |7|. 

Another important issue is that within MSSM inflation the quantum stability of 
the inflationary potential can be analyzed correctly as the couplings of the inflaton are 
known [|^. Since inflation is driven near the saddle point, it becomes an important 
question to ask how stable the potential is under large radiative corrections? To 
our surprise what we found is that an existence of a saddle point is not ruined 
by radiative corrections, although the point of inflection does shift towards higher 
vacuum expectation values (VEVs). The predictions for CMB also does not modify 
at all within the current uncertainties. The model predictions are also robust as there 
is no supergravity (SUGRA) corrections which can spoil the cosmological flatness of 
the potential. Although one would expect that SUGRA to play an important role, 
but all such corrections are absorbed when the choice of a saddle point is made 
where slow roll inflation occurs. In this respect one can take a view that a saddle 
point inflation can address the SUGRA-eta problem Of course SUSY is the 
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key ingredient which helps maintaining the flatness of the potential, in this respect 
SUSY along with gauge symmetry not only guides us towards a more holistic model 
of inflation, but also with sharp predictions which can be tested at the LHC and 
in future CMB experiments. Note that the model does not generate observationally 
significant stochastic gravitational wave background radiation or the cosmic strings. 

An interesting twist comes from string theory alone. As it is fairly well estab- 
lished by now that there exists a string theory landscape with plethora of vacua, for 
a review Therefore, in this landscape, only way our patch of the observable 

universe evolves to its present state regardless of how it began, if we only secure that 
the last stage of inflation were driven within an MSSM vacuum. In any case the 
MSSM fields (and couplings) must be in the low energy spectrum. At some stage 
the MSSM sector should take over. False vacuum inflation in MSSM then makes any 
previous stage in the history of the universe oblivious. However the observable part 
of the universe can still emerge since false vacuum inflation also sets the stage for a 
last bout of successful inflation (also happening in the MSSM sector). 

In the following sections we will discuss various virtues of MSSM inflation, in 
section 2, we discuss general principles behind a successful inflation, in section 3, we 
discuss MSSM and introduce a concept behind F-and D- flat directions. In section 
4, we describe MSSM inflation. In section 5, we discuss radiative and SUGRA 
corrections. In section 6, we discuss reheating and thermalization. In section 7, 
we discuss thermal neutralino cold dark matter creation. In section 8, our effort 
is to embed MSSM inflation with grand unified theories with an inclusion of the 
Majorana neutrinos. In section 9, we discuss how a small Yukawa coupling required to 
explain the neutrino masses could also maintain the flatness of the MSSM inflationary 
potential. In the same section we describe MSSM inflation within gauge mediated 
SUSY breaking scenarios. In the last section 10, we discuss initial conditions for an 
MSSM inflation, we speculate the role of string landscape in order to connect high 
and low scale inflations. We also point out a graceful exit from a string landscape 
which can happen via a slow roll MSSM inflation. 

2. What should be the inflation properties? 

In a minimal model, the inflaton is the only source of density perturbations and also 
responsible for a successful reheating. This is also true in the case of a curvaton [|60|] . 
Where the success of curvaton lies within a successful reheating of the SM degrees 
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of freedom during the curvaton dominance. It is possible to find an MSSM curvaton 
candidate which has all the properties to create an universe with the right relativistic 
degrees of freedom, cold dark matter and the seed for the primordial density pertur- 
bations pTI]. There are also variants of curvaton in the context of inhomogeneous 



reheating and generating density perturbations |6^. In either cases it is important 



to realize that a successful reheating ensures a successful CMB predictions, as an 
example in the case of an MSSM curvaton, recently we have singled out an unique 
MSSM flat direction which would create the right CMB predictions and the SM 
baryons along with neutralino type cold dark matter . 



The dream has been to embed the model of inflation in particle physics naturally. 
The first attempts were made by Guth [Q, where inflation would occur in a false 
vacuum of the Higgs of the Grand Unified Theory (GUT), fell into trouble, because 
the end of inflation happened due to first order phase transition. The bubble of true 
vacuum remains cold in the sea of false vacuum and there is no way to generate 
thermal entropy other than colliding the bubbles ( the potential energy is stored in 
the bubble walls ). There were attempts to address this issue, but all such attempts 
were not so attractive, as they would all involve fields such as absolute gauge singlets 



from a hidden sector, see for a review 12 



A successful reheating must ensure that: 

• SM Baryons are excited: 

The SM baryons are excited dominantly as light degrees of freedom to ensure a 
successful BBN . Note that BBN puts stringent constraints on hidden light 
degrees of freedom [?]. At best one could accommodate one light neutrino or 



any relativistic species . This suggests that the hidden sector particles must 
be heavier than few ~ TeV so that they kinetically and thermally decouple 
from the rest of the thermal bath before they could decay through Planck 
suppressed interactions. 

• Cold dark matter is created: 

A successful reheating must create conditions for generating thermally/non- 
thermally cold dark matter essential for the structure formation ||1CI|| . 

These two conditions do not necessarily bar the inflaton to be a part of a hidden 
sector. However if our main aim is to seek a model of particle physics then the 
hidden sector serves very little progress in our quest. Note that there is no underlying 
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symmetry which helps understanding the masses of the hidden sector fields and their 
couplings to the SM degrees of freedom. 

If any model of infiation seeks to be successful, then it must possess: 

• Credibility: 

The model parameters, such as mass, couplings are not chosen ad-hoc to match 
the CMB observations, rather they should arise naturally from an underlying 
theory. 

• Stability: 

As it is well known a slow roll infiation needs a fiat potential. By construction 
the infiaton energy density couples to gravity and the infiaton couples to at least 
the SM degrees of freedom, therefore, one needs to ensure that the background 
geometry, quantum corrections, supergravity effects do not spoil the flatness of 
the potential. 

• Testability: 

It would be desirable to test the micro physical ingredients of inflation in a 
terrestrial laboratory. It is fair to say that we have indeed tested several ideas 
of inflation in CMB physics, but in order to really seek the true origin of 
inflation one must do more than that. 

Based on the above issues, it is arguably simpler if inflation were driven solely 
by the SM particles. However SM spectra is full of fermions, only scalar is the SM 
Higgs but with a relatively small VEV to match the observational data 

This would eventually push the inflaton candidature to the physics beyond the 

SM. However this is the frontier where we lack our grounds on the experimental front, 

hence we are forced to the speculations. This is perhaps one of the main reasons 

for introducing absolute gauge singlets as the inflaton, or modeling the inflaton in a 

hidden sector. However, here as we have reiterated, we wish to have a concrete model 

of inflation which has a better predictions from the CMB observations, successful 

reheating and of course the model should be testable in future collider. 

^Could SM Higgs be the inflaton ? 
The Higgs searches at LEP has pushed the Higgs mass above 114 GeV |6^, this means that the 
electroweak phase transition via the SM Higgs is of second order or perhaps crossover in nature. 
Inflation could potentially work out if the Higgs field rolls very slowly, however, the Higgs VEV is 
sufficiently low enough to generate the scalar density perturbations to match the observed temper- 
ature anisotropy. 
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3. Supersymmetry as a tool 

Supersymmetry (SUSY) has many virtues, foremost, it is the most attractive sce- 
nario to address the hierarchy between the Planck and the electroweak scale. One 
particular advantage of SUSY is that the quantum corrections due to bosonic and 
fermionic loops cancel exactly in a SUSY limit, rendering the stability of masses, cou- 
plings and the scalar potentials. In order to address the hierarchy problem the SUSY 
must be broken spontaneously. At low energies, within MSSM, the SUSY breaking 
effects are captured by soft parameters, i.e. mass, trilinear couplings, etc. In low 
energy SUSY breaking scenarios, the scalar mass is ~ 0{TeV). Further note that 
SUSY allows many scalar fields (corresponding to every quarks and leptons within 
MSSM), therefore, it is interesting to ask what is the role of sfermions in cosmology? 

In softly broken SUSY, the quantum corrections give rise to Logarithmic run- 
ning to masses and couphngs. This ensures that the scalar potential is at least Log- 
arithmically stable under quantum corrections. In this respect the SUSY inspired 
inflationary potentials are at least stable under quantum corrections. This leads to 
satisfying one of the cornerstones of a successful inflationary model building. 

3.1 MSSM and its potential 

Let us remind the reader that the matter fields of MSSM are chiral superfields $ = 
(p + ^/29ip + 66F, which describe a scalar 0, a fermion ip and a scalar auxiliary field 
F. In addition to the usual quark and lepton superfields, MSSM has two Higgs fields, 
Hu and H^. Two Higgses are needed because W, which in the Standard Model gives 
masses to the w-quarks, is forbidden in the superpotential. 
The superpotential for the MSSM is given by |p6| 

Wmssm = XuQHuU + XdQHdd + XeLHae + ^iHuHd , (3.1) 

where Hu, Hd, Q, L, u, d, e in Eq. (|3.1|) are chiral superfields, and the dimensionless 
Yukawa couplings A^, A^, Ag are 3x3 matrices in the family space. We have suppressed 
the gauge and family indices. Unbarred fields are SU{2) doublets, barred fields 
SU{2) singlets. The last term is the /i term, which is a supersymmetric version of 
the SM Higgs boson mass. Terms proportional to H*Hu or H^Hd are forbidden in 
the superpotential, since Wmssm niust be analytic in the chiral fields. Hu and Hd 
are required not only because they give masses to all the quarks and leptons, but also 
for the cancellation of gauge anomalies. The Yukawa matrices determine the masses 
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and CKM mixing angles of the ordinary quarks and leptons through the neutral 
components of Hu = (-ff^,-ff°) and Ha = {H^,H^). Since the top quark, bottom 
quark and tau lepton are the heaviest fermions in the SM, we assume that only the 
(3, 3) element of the matrices A„, A^, Ag are important. In this limit only the third 
family and the Higgs fields contribute to the MSSM superpotential. 

The SUSY scalar potential V is the sum of the F- and D-terms and reads 

V = J2\^^\" + lj29lD^^^ (3.2) 

i a 

where 

p^^dWMSSM^ ]ja^^^j.a^_ (33) 

Here we have assumed that (pi transforms under a gauge group G with the generators 
of the Lie algebra given by T'*. 

3.2 F-and D-fiat directions 

For a general supersymmetric model with chiral superfields Xj, it is possible to 

find out the directions where the potential Eq. |3j^ vanishes identically by solving 
simultaneously 

X dW 

jya ^ ^^r^a^ ^ Q ^ Fx, = T;-r = 0- (3.4) 

0X4 



Field configurations obeying Eq. (|3.4|) are called respectively D-fiat and F-fiat. 

D-flat directions are parameterized by gauge invariant monomials of the chiral 
superfields. A powerful tool for finding the fiat directions has been developed in 
[p5| , |6^ , where the correspondence between gauge invariance and fiat directions has 
been employed. The configuration space of the scalar fields of the MSSM contains 
49 complex dimensions (18 for Qi, 9 each for Ui and di, 6 for Li, 3 for Cj, and 2 each 
for Hu and Hd), out of which there are 12 real D-term constraints (8 for SU{3)c, 3 
for SU{2)l, and 1 for f/(l)y), which leaves a total of 37 complex dimensions |6^, |6^] . 
The trick is to construct gauge invariant monomials forming SU{3)c singlets and 
then using them as building blocks to generate SU{3)c x SU{2)l, and subsequently 
the whole SU{3)c x SU{2)l x U{1)y invariant polynomials |6^, |66|. However these 
invariant monomials give only the D-fiat directions. For F-flat directions, one must 
solve explicitly the constraint equations Fx^ = 0. 

A single fiat direction necessarily carries a global U{1) quantum number, which 
corresponds to an invariance of the effective Lagrangian for the order parameter 
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B-L 


HuHd 





LHu 


-1 


udd 


-1 


QLd 


-1 


LLe 


-1 


QQud 





QQQL 





QLue 





uude 





QQQQu 


1 


QQuue 


1 


LLddd 


-3 


uuuee 


1 


QLQLdd 


-2 


QQLLdd 


-2 


uudddd 


-2 


QQQQdLL 


-1 


QLQLQLe 


-1 


QLuQQdd 


-1 


uuuddde 


-1 



Table 1: Renormalizable F and D flat directions in the MSSM 



under phase rotation (f) — > e*^0. In the MSSM the global U{1) symmetry is -B — L. 
For example, the Lif„-direction (see below) has B — L = —1. 

A flat direction can be represented by a composite gauge invariant operator, 
formed from the product of k chiral superfields $i making up the fiat direction: 
Xm = $i<I>2 • ■ ■ ^rn- The scalar component of the superfield Xm is related to the order 
parameter cj) through Xm = ccjf^. For a flat direction represented by polynomial the 



description is much more involved, see [37 



3.2.1 An example of F-and D-fiat direction 

The fiat directions in the MSSM are tabulated in Table 1. An example of a D-and 
F-fiat direction is provided by 




(3.5) 



where is a complex field parameterizing the fiat direction, or the order parameter, 
or the AD field. All the other fields are set to zero. In terms of the composite gauge 
invariant operators, we would write X^ = LHu {m = 2). 

From Eq. ( [4.14|) one clearly obtains F^^ = XuQu+fiHd = Fl = XdH^e = for all 



(j). However there exists a non-zero F-component given by F^^ = fiHu- Since fi can 
not be much larger than the electroweak; scale M\y ~ (9(1) TeV, this contribution 
is of the same order as the soft supersymmetry breaking masses, which are going to 
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lift the degeneracy. Therefore, following [55|, one may nevertheless consider LHu to 
correspond to a F-flat direction. 
The relevant D-terms read 

D%um = Hir.H^ + LtraL = _ ^ o . (3.6) 

Therefore the LH^ direction is also D-flat. 

The only other direction involving the Higgs fields and thus soft terms of the 
order of /i is HuH^. The rest are purely leptonic, such as LLe, or baryonic, such as 
udd, or mixtures of leptons and baryons, such as QLd. 

3.2.2 Lifting by non-renormalizable operators 

Non-renormalizable superpotential terms in the MSSM can be viewed as effective 
terms that arise after one integrates out fields with very large mass scales appearing 
in a more fundamental (say, string) theory. Here we do not concern ourselves with 
the possible restrictions on the effective terms due to discrete symmetries present in 
the fundamental theory, but assume that all operators consistent with symmetries 
may arise. Thus in terms of the invariant operators Xm, one can have terms of the 



type 38, 35 



where the dimensionality of the effective scalar operator d = mk, and h is a coupling 
constant which could be complex with \h\ ~ Here M is some large mass, 

typically of the order of the Planck mass or the string scale (in the heterotic case 
M ~ Mqut)- The lowest value of is 1 or 2, depending on whether the fiat direction 
is even or odd under i?-parity. 

A second type of term lifting the fiat direction would be of the form [pSl 



where ip is not contained in X^- The superpotential term Eq. ( |3.8| ) spoils F- flatness 
through ^ 0. An example is provided by the direction u\U2Uze\e2-, which is 
lifted by the non-renormalizable term W = {h' / M)uiU2d2ei. This superpotential 
term gives a non-zero contribution = [h' / M)uiU2ei ~ {h'/M)(j)^ along the fiat 
direction. 

Assuming minimal kinetic terms, both types discussed above in Eqs. (|3.7]j3.^ ) 
yield a generic non-renormalizable potential contribution that can be written as 



lA 



2 
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Figure 1: The colored curves depict the full potential, where V{x) = 
V{(j))/{0.5 m^M|(m^/Mp)i/2)^ ^ = (A„Mp/m0)i/'^(0/Mp). The black curve is the 
potential arising from the soft SUSY breaking mass term. The black dots on the colored 
potentials illustrate the gradual transition from minimum to the saddle point and to the 
maximum. 



where we have defined the coupling |Ap = + By virtue of an accidental 

i?-symmetry under which has a charge R = 2/d, the potential Eq. (p.9|) conserves 
the U{1) symmetry carried by the fiat direction, in spite of the fact that at the 
superpotential level it is violated, see Eqs. ( |3.7| , p.8| ). The symmetry can be violated if 
there are multiple fiat directions, or by higher order operator contributions. However 
it turns out that the B — L violating terms are always subdominant. This is 
of importance for baryogenesis considerations, where the necessary B — L violation 
should therefore arise from other sources, e.g. such as soft supersymmetry breaking 
terms. 



4. Gauge Invariant Inflaton 



Let us recapitulate the main features of MSSM flat direction inflation [Q, ^ 
The framework is solely based on MSSM together with gravity, so consistency dictates 
that all non-renormalizable terms allowed by gauge symmetry and supersymmetry 
should be included below the cut-off scale, which we take to be the Planck scale. 
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The superpotential term which hfts the F-flatness is given by: 

n>3 

where $ is a gauge invariant superfield which contains the flat direction. Within 
MSSM all the flat directions are lifted by non-renormalizable operators with 4 < 
n < 9 where n depends on the flat direction. We expect that quantum gravity 
effects yield M = Mp = 2.4 x 10^^ GeV and A„ ~ C(l) Note however that 

our results will be valid for any values of A„, because rescaling A„ simply shifts the 
VEV of the flat direction. Let us focus on the lowest order superpotential term in 
Eq. ( |9.1| ) which lifts the flat direction. Soft SUSY breaking induces a mass term for 
and an A-term so that the scalar potential along the flat direction reads 

V = -™j + A coB(n9 + + , (4.2) 

Here and 9 denote respectively the radial and the angular coordinates of the com- 
plex scalar field $ = exp[i^], while 9 a is the phase of the A-term (thus A is a 
positive quantity with dimension of mass). Note that the first and third terms in 



Eq. (|9.17|) are positive definite, while the A-term leads to a negative contribution 
along the directions whenever cos{n9 + 9a) < 0, see |l| ^. There are other attempts 
to embed inflation within a gauge theory they however do not explain SM 
reheating, besides they also have some caveats related to Planckian VEVs. 

4.1 Inflation near the Saddle Point 

The maximum impact from the A-term is obtained when cos{n9 + 9a) = — 1 (which 
occurs for n values of 9). Along these directions V has a secondary minimum at 
= 00 ~ {'m<j,Mp^'^Y^"' ^ -C Mp (the global minimum is at = 0), provided that 

A2>8(n-l)mJ. (4.3) 

At this minimum the curvature of the potential is positive both along the radial 
and angular directions. If the A is too large, the secondary minimum will be deeper 
than the one in the origin, and hence becomes the true minimum. However, this 



^The importance of the A-term was first highlighted in a successful MSSM curvaton model |63| , 
where again the curvaton carries the SM charges. By implementing so it also leads to a successful 
CMB prediction, SM reheating and a detectable signature at the LHC. 
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is phenomenologically unacceptable as such a minimum will break charge and/or 
color HI, H. With a total potential: V ~ m^^g ^ ^2 (rri^M^-^)^^^""^^ 

As discussed in if the local minimum is too steep, the field will become 

trapped there with an ensuing infiation that has no graceful exit like in the old infia- 
tion scenario On the other hand in an opposite limit, with a point of infiection, 
a single fiat direction cannot support infiation [^, one would require an assisted 



infiation with the help of many fiat directions fl^ .. 

However, in the gravity mediated SUSY breaking case, the A-term and the soft 
SUSY breaking mass terms are expected to be the same order of magnitude as the 
gravitino mass, i.e. 



m,/) ~ A ~ m3/2 ~ C>{1) TeV 



(4.4) 



Therefore, as pointed out in [Q, in the gravity mediated SUSY breaking it is pos- 
sible that the potential barrier actually disappears and the inequality in Eq. ( |9.5|) is 
saturated so that A and are related by 



^2 = - l)ml 



(4.5) 



If the above condition is satisfied then both the first and second derivatives of V 
vanish at 0O; i-e. V'^cpo) = 0, V'(0o) = 0. As the result, if initially ~ 0o; a slow 
roll phase of infiation is driven by the third derivative of the potential. 

Note that this behavior does not seem possible for other SUSY breaking scenarios 
such as the gauge mediated breaking or split SUSY [^]. In split SUSY the A- 
term is protected by an i?-symmetry, which also keeps the gauginos light while the 
sfermions are quite heavy []72| ^. 

4.2 Slow roll 

The potential near the saddle point Eq. ( ^.5| ) is very fiat along the real direction but 
not along the imaginary direction. Along the imaginary direction the curvature is 
determined by m^. Around (po the field lies in a plateau with a potential energy 

2)2 



^(0o) 



[n 



2n{n 



2 ±2 



(4.6) 



^In the gauge mediated case there is an inherent mismatch between A and m^, except at very 



large field values where Eq. (4.4) can be satisfied. However there exists an unique possibility of a 
saddle point inflation within gauge mediated case which we will discuss in Section 9.2 |47[|. 
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with 



Xn\/2n - 2, 

This results in Hubble expansion rate during inflation which is given by 



_ (n-2) „.^^, 

v/6n(n - 1) Mp " ^^'"^ 

When (f) is very close to <^0; the flrst derivative is extremely small. The fleld is 
effectively in a de Sitter background, and we are in self-reproduction (or eternal 
inflation) regime where the two point correlation function for the flat direction fluc- 
tuation grows with time. But eventually classical friction wins and slow roll begins 

at ^ 0,elf H, H 

(0O-0self)-(^)'V (4.9) 

The regime of eternal inflation plays an important role in addressing the initial 
condition problem, see section 10. 

The observationally relevant perturbations are generated when 0cobe- 
The number of e-foldings between ^cobe and ^end, denoted by A/cobe follows from 

Eq. (IToleD 

^COBE ^ TTT^ T. (4.10) 

2n{n - l)M^(0o - 0cobe) 



The amplitude of perturbations thus produced is given by 



^--H^^-B^Vr'-^X-^K^)-^-. (4.11, 



and the spectral tilt of the power spectrum and its running are found to be P3, H5 



= l + 2r/-6e ~ 1 - -r^ — , (4.12) 

A/cOBE 

^^•^ ^ .. (4.13) 



dXnk A^cobe' 
4.3 Infiaton Properties and predictions 

4.3.1 Infiaton candidates 

As discussed in among nearly 300 flat directions there are two that can lead 

to a successful inflation along the lines discussed above. 

One is udd which, up to an overall phase factor, is parameterized by 

< = ^0, 4 = ^0, d^k = ^4>- (4.14) 
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Here 1 < a,/3,7 < 3 are color indices, and 1 < i, j, k < 3 denote the quark families. 
The flatness constraints require that a 7^ /? 7^ 7 and j 7^ k. 

The other direction is LLe ^, parameterized by (again up to an overall phase 
factor) 

where 1 < a, 6 < 2 are the weak isospin indices and 1 < i, j, k < 3 denote the lepton 
families. The flatness constraints require that a ^ h and i ^ j ^ k. Both these flat 
directions are lifted by = 6 non-renormalizable operators [|6^] , 

We 3 -^{LLe){LLe) , D ^{udd){udd) . (4.16) 

The reason for choosing either of these two flat directions* is twofold: (i) a non- 
trivial y4-term arises, at the lowest order, only at n = 6; and (ii) we wish to obtain 
the correct COBE normalization of the CMB spectrum. 

Those MSSM flat directions which are lifted by operators with dimension n = 7, 9 
are such that the superpotential term contains at least two monomials, i.e. is of the 
type 

W ^^<l>"-^ (4.17) 

If represents the flat direction, then its VEV induces a large effective mass term for 
ip, through Yukawa couplings, so that (ip) = 0. Hence Eq. ([4.17|) does not contribute 
to the y4-term. 

More importantly, as we will see, all other flat directions except those lifted by 
n = 6 fail to yield the right amplitude for the density perturbations. Indeed, as can 
be seen in Eq. ( |9.(j| ), the value of 0o, and hence also the energy density, depend on n. 

4.3.2 Infiaton Predictions 

According to the arguments presented above, successful MSSM flat direction inflation 
has the following model parameters: 

m<^~l-10TeV, n = 6, A = V^m^ , A ~ 0(1) . (4.18) 
^When the flat direction develops a VEV during inflation, it spontaneously breaks SU (2) x U{l)y, 

which gives masses to the corresponding gauge bosons. It is possible to obtain a seed perturbations 

for the primordial magnetic field in this case, see [Q. 

^ Since LLe are uddq are independently D- and i^-flat, inflation could take place along any of 

them but also, at least in principle, simultaneously. The dynamics of multiple flat directions are 

however quite involved [p7|. 
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Here we assume that A (we drop the subscript "6") is of order one, which is the most 
natural assumption when M = Mp. 

The Hubble expansion rate during inflation and the VEV of the saddle point 

9 

are 

i/inf -1-10 GeV , 00 ~ (1 - 3) X 10^^ GeV . (4.19) 

Note that both the scales are sub-Planckian. The total energy density stored in the 
inflaton potential is Vq ~ 10^^ — 10'^^ GeV^. The fact that 0o is sub-Planckian guar- 
antees that the inflationary potential is free from the uncertainties about physics at 
super-Planckian VEVs. The total number of e-foldings during the slow roll evolution 



is large enough to dilute any dangerous relic away [45 



Kot ~ 10' , (4.20) 

At such low scales as in MSSM inflation the number of e-foldings, A/cobe, required 
for the observationally relevant perturbations, is much less than 60 |71| , [75| . If the 
inflaton decays immediately after the end of inflation, we obtain A/cobe ~ 50. De- 
spite the low scale, the flat direction can generate adequate density perturbations 
as required to explain the COBE normalization. This is due to the extreme flatness 
of the potential (recall that V = 0), which causes the velocity of the rolling flat 



direction to be extremely small. From Eq. (^4.11|) we flnd an amplitude of 



<5h ^ 1-91 X 10-^ (4.21) 

There is a constraint on the mass of the flat direction from the amplitude of the 
CMB anisotropy: 

^ (100 GeV) x (^^%^) ' • (4-22) 

We get a lower limit on the mass parameter when A < 1. For smaller values of 
A ^ 1, the mass of the flat direction must be larger. Note that the above bound on 
the inflaton mass arises at high scales, i.e. (p = 0o- However, through renormalization 
group flow, it is connected to the low scale mass, as will be discussed in Sect. 4. 

The spectral tilt of the power spectrum is not negligible because, although the 
flrst slow roll parameter is e ~ I/A/'cqbe ^ l' other slow roll parameter is given 



^We note that i?inf and 0o depend very mildly on A as they are both cx A 
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by ^7 = — 2/A/coBE and thus, see Eq. ( |4.12|) 



10 



- 0.92, (4.23) 
-^--0.002, (4.24) 

where we have taken A/cobe ~ 50 (which is the maximum value allowed for the scale 
of inflation in our model). In the absence of tensor modes, this agrees with the current 



WMAP 3-years' data within 2a |T0|. Note that MSSM inflation does not produce 
any large stochastic gravitational wave background during inflation. Gravity waves 
depend on the Hubble expansion rate, and in our case the energy density stored in 
MSSM inflation is very small. 

4.4 Departure from the saddle point 

Inflation can still happen for small deviations from the saddle point condition Eq. ( [4.51 ). 
To quantify this, we define a parameter such that [^, Q : 

a\ (4.25) 



8{n - \)m\ V 2 

For o? 7^ 0, the saddle point becomes a point of inflection where V^"(0o) = 0, and 

"a - 2\2 



^'(0o) = \-^) o^^i^o. (4.26) 

If < 0, the potential has a local minimum and a maximum. In this case the 
flat direction is trapped in the local minimum. It will eventually tunnel past the 
maximum and a period of slow roll inflation will follow [|^. If > 0, the potential 
has no maximum or local minimum, and then slow roll inflation occurs around 0o- 
For ^ the expressions for and 6h are modified as (see also |5 



1/2 .m^Mp 1 . 2r 



and 



where 



^ 5^V 3''^'' " ^^^"^ " 0r sin^iATcoBEVA^] , (4.27) 



= 1 - 4v^ cotiATcoBEV^], , (4.28) 



A2^n2(n-l)VAr^oBE(^)'- (4.29) 



10 



Obtaining > 0.92 (or < 0.92, which is however outside the 2a aUowed region) requires 



deviation from the saddle point condition in Eq. (4.5), see the next subsection. For a more detailed 



discussion on the spectral tilt, see also Rcfs. [[76|,[p0[ 
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Figure 2: is plotted as a function of for different values of m^. A is defined in the 
text. We choose A =1. 



Note that for for = 0, Eqs. (|427|JOg ) are reduced to (|4Tl|j4rT^ ) respectively. 
For < 0, the spectral index will be smaller than that in Eq. ( [4.12| ), thus outside 
the 2a region from observations. The more interesting pointed out in [BO], 

happens for > 0. We can in this case get all values within the allowed range 



0.92 <n,<l 111 for 



< < 



TT 



4A/coBE 



(4.30) 



The inflaton mass, m^, is constrained by the experimental data on the spectral 
index rig [T^, [Tl| and 6h fT^ • 

We first find the solutions of by solving Eqs. ([4.27| , [4.28|) . Ug depends mainly 
on and is mostly independent of and A (the coupling in Eq. (1)). The param- 
eter A^ is defined in Eq. ( [4.30| ). We therefore solve A^ from Eq. ( |4.28| ) and apply 
this solution to determine the bounds on from the Eq. (^4.27| ). In figure 0, we 
show function of A^. The range for A^ is determined from Eq. ( |4.30| ). 



In figure ^, we show Sh as a function of for different values of m^. The blue 
band shows the experimentally allowed region. We find that smaller values of 
are preferred for smaller values of n^. We also find that the allowed range of is 
75 — 440 GeV for the experimental ranges of and 6h- We assume A ~ 1 for these 
two figures. If A is less than 0{1), e.g., A ~ 0.1 or so (which can occur in 5*0(10) 
model), it will lead to an increase in m^. We will need to study these allowed ranges 
of the inflaton mass in the mSUGRA model. Since the inflaton mass is related to the 



- 22 - 



xlO 



0.3 



CO 



0.25 



^ 0.2 



0.15 



0.1 



I I I I I I I I I U I I I 



I I I I I 



0.925 0.95 0.975 1 

n 



Figure 3: 6h is plotted as a function of for different values of m^. We used A =1. The 
blue band denotes the experimentally allowed values of 6h- 

parameters of the mSUGRA model, the main question is whether the allowed range 
of the inflaton mass is consistent with the experimentally allowed mSUGRA model 
or not. See section 7. 



5. Radiative and supergravity corrections 

Since the MSSM inflaton candidates are represented by gauge invariant combina- 
tions which are not singlets. The inflaton parameters receive corrections from gauge 
interactions which, unlike in models with a gauge singlet inflaton, can be computed 
in a straightforward way. Quantum corrections result in a logarithmic running of the 
soft supersymmetry breaking parameters and A. In this section we will discuss 
running of the potential with VEV-dependent values of m^^cj)) and in Eq. (^4.5| ). 
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Our conclusion is that the running of the gauge couphngs do not spoil the existence of 
a saddle point. However the VEV of the saddle point is now displaced; by how much 
will depend precisely on the inflaton candidate. In order to discuss the situation, we 
derive a general expression for the one-loop effective potential for the flat directions, 
and then focus on the LLe direction, for which the system of Renormalization Group 
(RG) equations can be solved analytically 



5.1 One-loop effective potential 



The first thing to check is whether the radiative corrections remove the saddle point 
altogether. The object of interest is the effective potential at the phase minimum 



^^min = TT, for which we obtain [H5 



2 J 2 

-mo0 



1 + i^l log J 2 



X2 



/>2(n-l) 



1 + l0£ 



2' 



1 + K2 log 



(5.1) 



where mo, Aq, and A„,o are the values of m^, A and A„ given at a scale yUo- Here Aq 
is chosen to be real and positive (this can always be done by re-parameterizing the 
phase of the complex scalar field 0), and \Ki\ < 1 are coefficients determined by the 
one-loop renormalization group equations. 

Our aim is to find a saddle point of this effective potential, so we calculate the 1st 
and 2nd derivatives of the potential and set them to zero. This is a straightforward 



although somewhat cumbersome exercise that results in the expression 45 



n-3 



4A„(n - 1 + 



A[l+ '^K^ ± W A^ (l + ^K2] - 8m2(l + K^){n-1 + Ks) 



(5.2) 

where m^. A, and A„ are values of the parameters at the scale (po- Inserting this into 



^udd case follows similar discussion but requires numerics to solve the equations. 
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V(f)<t> = 0, we can then find the condition to have a saddle point at (j)o: 



= 2ml(n - 1 + K^)F^F2F^ 



F2 
F3 



l + Ki 



(n - l){2n - 3) + {An - 5)K3^ 
-(l + SiTi) 



n ^ 

n-l + Ks 



(n - l)(2n - 3) + {An - 5)K3^ - (n - 1 + 



1 + -K, 

n 



n 



(5.3) 



In the limit when \Ki\ <^ 1, this mercifully simplifies to 45 



= 8(n - l)m^ 



M^-^m^{(j)o) 



1 + K^--K2 + 

n n — 1 



K^2{n - 1) 



2 



1 



2(n - 1 



(5.4) 
(5.5) 



Note that Eqs. (O.^) give the necessary relations between the values of and 
A as calculated at the saddle point. The coefficients Ki need to be solved from the 
renormalization group equations at the scale given by the saddle point /i = 0o- Since 
Ki are already one loop corrections, taking the tree-level saddle point value as the 
renormalization scale is sufficient. 

The conclusion is robust, although the soft terms and the value of the saddle 
point are all affected by radiative corrections, they do not remove the saddle point 
nor shift it to unreasonable values. The existence of a saddle point is thus insensitive 
to radiative corrections. 



5.2 RG equations for the LLe direction 

The form of the relevant RG equations depend on the fiat direction. RG equations 
for LLe are simpler since only the SU{2)w x U{1)y gauge interactions are involved 
and the lepton Yukawa couplings are negligible. The case of udd requires numerics 
if u is chosen from the third family. For other choices, however, it closely resembles 
LLe. For LLe the one-loop RG equations governing the running of m|. A, and A 
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with the scale /i are given by 



u , = TT -rrirtOr, H — mififi 

ciA 1 /3 . 2 3 . 2 



d/i 27r2 V2 2 



Here mi, m2 denote the mass of the U{1)y and SU{2)\y gauginos respectively and 
gi, g2 are the associated gauge couplings. It is a straightforward exercise to obtain 
the equations that govern the running of A and A associated with the [LLe)'^ su- 
perpotential term (which lifts the LLe fiat direction). Note that L has the same 
quantum numbers as H^, and hence in this respect LLe combination behaves just 
like HdLe. One can then use the familiar RG equations that govern the Yukawa 
coupling and A-term associated with the H^Le superpotential term [Q. However, 



as explained in |7^, the coefficients of the terms on the right-hand side are propor- 
tional to the number of superfields contained in a superpotential term. Hence the 
second and third equations in (|5.6| ) are simply obtained from those for the H^Le term 
after multiplying by a factor of 2. The first equation in ( |5.6| ) is also easily found by 
taking the electroweak charges of Lj, Lj and e superfields into account while taking 
into account that m| = (m|. + m|. + m^)/3. 



The running of gauge couplings and gaugino masses obey the usual equations [pq] - 



dgi 


11 3 


dfi 




dg2 


1 3 


dfi 




rhi\ 


d fm2 


9l) 


df^ V 9i 



^l^|-.^l^|-0. (5.7) 



^We would like to thank Manuel Drees for explaining this point to us. 
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The solutions of the RG equations are |45 

9iM 



9i 



rrii 



1 



A 
A 



v4(/io) + 6 (m2(/io) - ^2) + — (mi(/io) - mi) 



A(/io 



92 



9i 



(5.8) 

(5.9) 

(5.10) 
(5.11) 

(5.12) 



where i = 1,2, 61 = II/Stt^ and 62 = I/Stt^. Ignoring the running of the gaugino 



masses and gauge couphngs, we find that 

1 



47r2 
3 

3 

'8^ 



m2 



9I + 



mi 



9l 



Ao [9I + 9f\ 



(5.13) 



where the subscript denotes the values of parameters at the high scale /iq. 

For universal boundary conditions, as in minimal grand unified supergravity, 
the high scale is the GUT scale fix ~ 3 x 10^^ GeV, rhi{^x) = 'fn2{fJ'x) = m 
and gi = a/vt/IO ~ 0.56, g2 = a/vt/G ~ 0.72. Then we just use RG equations to 
run the coupling constants and masses to the scale of the saddle point /io = 0o ~ 
2.6 X 10^^ GeV for Mp = 2.4 x 10^^ GeV, m^, = 1 TeV, Aq = 1. With these values 
we obtain 



K2 
K3 



-o.oi7r, 

-0.0085^, 
-0.029. 



(5.14) 
(5.15) 
(5.16) 



where ^ = fh/m^ is calculated at the GUT scale. 

Typically the running based on gaugino loops alone results in negative values of 
Ki [^. Positive values can be obtained when one includes the Yukawa couplings, 
practically the top Yukawa, but the order of magnitude remains the same. 

Thus radiative corrections modify a and we need to fine tune the potential to a 
few (but not all) orders in perturbation theory. 
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Figure 4: The running of for the LLe inflaton when the saddle point is at (f)Q = 
2.6 X lO^^GeV (corresponding to n = 6, = 1 TeV and A = 1). The three curves 
correspond to different values of the ratio of gaugino mass to flat direction mass at the 
GUT scale: ^ = 2 (dashed), ^ = 1 (sohd) and i = 0.5 (dash-dot). 

5.3 Aq vs. 

One final comment is in order before closing this Section. Unlike m^, there is no 
prospect of measuring the A term, because it is related to the non-renormalizable 
interactions which are suppressed by Mp. However, a knowledge of super symmetry 
breaking sector and its communication with the observable sector may help to link 
the non-renormalizable A-term under consideration to the renormalizable ones. 

To elucidate this, let us consider the Polonyi model where a general A-term at 
a tree level is given by 

m^,2[{a-?,)W + <P{dW/d(t))l 

with a = 3 — -\/3 . One then finds a relationship between A-terms corresponding 
to n = 6 and n = 3 superpotential terms, denoted by Aq and ^3 respectively, at high 
scales: 

A, = l^A,. (5.17) 
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One can then use relevant RG equations to relate Aq which is relevant for inflation, 
to A3 at the weak scale, which can be constrained and/or measured. In principle 
this can also be done in general, provided that we have sufficient information about 
the supersymmetry breaking sector and its communication with the MSSM sector. 



see some related discussions in 81 



5.4 Supergravity corrections 

SUGRA corrections often destroy the slow roll predictions of inflationary potentials; 
this is the notorious SUGRA-?] problem [^. In general, the effective potential 



depends on the Kahler potential K as V ^ ^e^'-'^*''^''/*^pV^(0) j so that there is a 
generic SUGRA contribution to the fiat direction potential of the type 



= H'M^f j , (5.18) 

where / is some function (typically a polynomial). Such a contribution usually gives 
rise to a Hubble induced correction to the mass of the fiat direction with an unknown 
coefficient, which depends on the nature of the Kahler potential 

Let us compare the non-gravitational contribution, Eq. ( p.l7| ), to that of Hubble 
induced contribution, Eq. ( [5.181 ). Writing / ~ (0/Mp)^ where p > 1 is some power. 



we see that non-gravitational part dominates whenever 

^mfM| (^^y « m|02 , (5.19) 

so that the SUGRA corrections are negligible as long as </)o <C Mp, as is the case 
here (note that ifinfMp ~ m^cf)^). The absence of SUGRA corrections is a generic 
property of this model. Note also that although non-trivial Kahler potentials give 
rise to non-canonical kinetic terms of squarks and sleptons, it is a trivial exercise to 
show that at sufficiently low scales, ifjnf ^ m^, and small VEVs, they can be rotated 
to a canonical form without affecting the potential 



^^If the Kahler potential has a shift symmetry, then at tree level there is no Hubble induced correc- 
tion. However, at one-loop level relatively small Hubble induced corrections can be induced js^, |83|. 

^^The same reason, i.e. iJjnf ^ also precludes any large Trans-Planckian correction. Any 
such correction would generically go as (iJinf/A/»)^ ^ 1, where is the scale at which one would 
expect Trans-Planckian effects to kick in. 
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6. Reheating and Thermalization 



After the end of inflation, the flat direction starts rolling towards its global minimum. 
At this stage the dominant term in the scalar potential will be: m(^(/)^/2. Since the 
frequency of oscillations is ~ ~ lO'^ifinf, the flat direction oscillates a large 
number of times within the first Hubble time after the end of inflation. Hence the 
effect of expansion is negligible. 

We recall that the curvature of the potential along the angular direction is much 
larger than -ff^f. Therefore, the flat direction has settled at one of the minima along 
the angular direction during inflation from which it cannot be displaced by quantum 
fluctuations. This implies that no torque will be exerted, and hence the flat direction 
motion will be one dimensional, i.e. along the radial direction. 

Flat direction oscillations excite those MSSM degrees of freedom which are cou- 
pled to it. The inflaton, either LLe or udd flat direction, is a linear combination of 
slepton or squark fields. Therefore inflaton has gauge couplings to the gauge/gaugino 
fields and Yukawa couplings to the Higgs/Higgsino fields. As we will see particles with 
a larger couplings are produced more copiously during inflaton oscillations. There- 
fore we focus on the production of gauge fields and gauginos. Keep in mind that 
the VEV of the MSSM flat direction breaks the gauge symmetry spontaneously, for 
instance udd breaks SU{?))c x U{1)y while LLe breaks SU{2)v/ x t/(l)y, therefore, 
induces a supersymmetry conserving mass ~ g{(f){t)) to the gauge/gaugino fields in 
a similar way as the Higgs mechanism, where g is a gauge coupling. When the flat 
direction goes to its minimum, {(f){t)) = 0, the gauge symmetry is restored. In this 
respect the origin is a point of enhanced symmetry . 



There can be various phases of particle creation in this model, here we briefly 
summarize them below. Let us elucidate the physics, by considering the case when 
LLe flat direction is the inflaton 

• Tachyonic preheating: 

Right after the end of inflation, when we are close to the saddle point, the 
second derivative is negative. One might suspect that this would trigger tachy- 
onic instability in the inflaton fluctuations which will then excite the inflaton 
couplings to matter [BR . 
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^Reheating happens quickly due to a flat direction motion which is strictly one dimensional in 
our case. Our case is really exceptional, usually, the flat direction motion is restricted to a plane, 
which precludes preheating all together, for instance see pB, Bq, pTj. 
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However the situation is different in our case. As mentioned, only infiaton 
fluctuations witli a pliysical momentum k < will have a tachyonic instabil- 
ity. Moreover V" < only at field values which are ~ 0o- Tachyonic effects 
are therefore expected be negligible since, unlike the case in [^, the homo- 
geneous mode has a VEV which is hierarchically larger than (we remind 
that 00 > 10^'' GeV) and oscillates at a frequency u ~ m^. Further note 
fields which are coupled to the infiaton acquire a very large mass ~ h(f)Q from 
the homogeneous piece which suppresses non-perturbative production of their 
quanta at large infiaton VEVs. We conclude that tachyonic effects, although 
genuinely present, do not lead to significant particle production in our case. 

Instant preheating: 

An efficient bout of particle creation occurs when the infiaton crosses the origin, 
which happens twice in every oscillation. The reason is that fields which are 
coupled to the infiaton are massless near the point of enhanced symmetry. 
Mainly electroweak gauge fields and gauginos are then created as they have 
the largest coupling to the fiat direction. The production takes place in a short 
interval. At ~ {gm^(f)Q)~^^'^ , where 0o ~ 10^^ GeV is the initial amplitude 
of the infiaton oscillation, during which quanta with a physical momentum 

1 /2 

k <, ((^m^^o) are produced. The number density of gauge/gaugino degrees 
of freedom is given by |0D|], see also PT 



As the infiaton VEV is rolling back to its maximum value (po, the mass of the 
produced quanta g{(j){t)) increases. The gauge and gaugino fields can (pertur- 
batively) decay to the fields which are not coupled to the infiaton, for instance 
to (s)quarks. Note that (s)quarks are not coupled to the flat direction, hence 
they remain massless throughout the oscillations. The total decay rate of the 
gauge/gaugino fields is then given by F = C {g^ /A8tt) g(f), where C ~ (9(10) is 
a numerical factor counting for the multiplicity of final states. 

The decay of the gauge/gauginos become efficient when |^ 

(6.2) 

Here we have used (0(t)) ~ (j)om^t, which is valid when m<^t ^ 1, and F ~ 
t~^, where t represents the time that has elapsed from the moment that the 
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inflaton crossed the origin. Note that the decay is very quick compared with 
the frequency of inflaton oscillations, i.e. F ^ m^. It produces relativistic 
(s) quarks with an energy P5|] : 



E = -gm ^ -—^ . (6.3) 



2^"^ ^ V Cg 

The ratio of energy density in relativistic particles thus produced p^e^with re- 
spect to the total energy density po follows from Eqs. (|6.1| , |6.3|) : 



^-10-2^, (6.4) 
Po 

where we have used C ~ (9(10). This implies that a fraction ~ (9(10~^) of 
the inflaton energy density is transferred into relativistic (s) quarks every time 
that the inflaton passes through the origin. This is so-called instant preheating 
mechanism It is quite an efficient mechanism in our model as it can 



convert almost all of the energy density in the inflaton into radiation within a 
Hubble time (note that H^l ~ lO^m^^). 

6.1 Towards thermal equilibrium 

A full thermal equilibrium is reached when a) kinetic and b) chemical equilibrium 
are established. The maximum (hypothetical) temperature attained by the plasma 
would be given by: 

Tma. ~ V'/^ ~ (m^0o)'/' > 10^ GeV . (6.5) 

This temperature may be too high and could lead to thermal overproduction of grav- 
itinos However the dominant source of gravitino production in a thermal 



bath is scattering which include an on-shell gluon or gluino leg. In the next subsec- 
tion we describe a natural solution to this problem and show that the flnal reheat 
temperature is actually well below Eq. (|9.37|) , i.e. Tr <^ T^ax- 

One comment is in order before closing this subsection. The gravitinos can 
also be created non-perturbatively during inflaton oscillations, both of the helicity 
±3/2 and helicity ±1/2 states In models of high scale inflation (i.e. i^inf ^ 
m3/2) helicity ±1/2 states can be produced very efficiently (and much more copiously 
than helicity ±3/2 states). At the time of production these states mainly consist of 



-'^^In a favorable condition the fiat direction VEV coupled very weakly to the flat direction inflaton 
could also enhance the perturbative decay rate of the inflaton [BSl . 
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the inflatino (inflaton's superpartner). However these fermions also decay in the form 
of inflatino, which is coupled to matter with a strength which is equal to that of the 
inflaton. Therefore, they inevitably decay at a similar rate as that of inflaton, and 



hence pose no threat to primordial nucleosynthesis flOO . 

In the present case ~ ^7^3/2 ^ -f^inf- Therefore low energy supersymmetry 
breaking is dominant during inflation, and hence helicity ±1/2 states of the gravitino 
are not related to the inflatino (which is a linear combination of leptons or quarks) at 
any moment of time. As a result helicity ±1/2 and ±3/2 states are excited equally, 
and their abundances are suppressed due to kinematical phase factor. Moreover there 
will be no dangerous gravitino production from perturbative decay of the inflaton 



quanta ||9^, ^ |101|| . The reason is that the inflaton is not a gauge singlet and 



has gauge strength couplings to other MSSM fields. This makes the inflaton 
inflatino ± gravitino decay mode totally irrelevant. 

6.2 Solution to the gravitino problem 

In order to suppress thermal gravitino production it is sufficient to make gluon and 
gluino fields heavy enough such that they are not kinematically accessible to the 
reheated plasma, even if other degrees of freedom reach full equilibrium (for a detailed 



discussion on thermalization in supersymmetric models and its implications, see |^ 
|8^). This suggests a natural solution to the thermal gravitino problem in the case 
of our model. Consider another fiat direction with a non-zero VEV, denoted by y?, 
which spontaneously breaks the S'f/(3)c group. For example, if LLe is the infiaton, 
then udd provides a unique candidate which can simultaneously develop VEV ^^The 
induced mass for gluon/gluino fields will be: 

rriG ^ g{v{t)) < g(t)Q. (6.6) 

The inequality arises due to the fact that the VEV of if cannot exceed that of the 
inflaton cj) since its energy density should be subdominant to the inflaton energy 
density. 

If gifQ ^ Tjnax the gluon/gluino flelds will be too heavy and not kinematically 
accessible to the reheated plasma. Here ipo is the VEV of udd at the beginning 
of inflaton oscillations. In a radiation-dominated Universe the Hubble expansion 



"'^''To develop and maintain such a large VEV, it is not necessary that udd potential has a saddle 
point as well. It can be trapped in a false minimum during inflation, which will then be lifted by 
thermal corrections when the inflaton decays (as discussed in the previous subsection) 
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redshifts the flat direction VEV as oc H^^^^ which is a faster rate than the 
change in the temperature T oc H^^^ . Once g{ip) ~ T, gluon/gluino fields come into 
equilibrium with the thermal bath. As pointed out in Refs. |9^, ^ ^ , if the initial 
VEV of udd is 

(/?o>10^°GeV, (6.7) 

then the temperature at which gluon/gluino become kinematically accessible, i.e. 
g{Lp) ~ T, is given by [||] 

TR<10^GeV. (6.8) 

This is the final reheat temperature at which gluons and gluinos are all in thermal 
equilibrium with the other degrees of freedom. The standard calculation of thermal 
gravitino production via scatterings can then be used for T < Tr. Note however that 
Tr is sufficiently low to avoid thermal overproduction of gravitinos. 

Finally, we also make a comment on the cosmological moduli problem. The 
moduli are generically displaced from their true minimum if their mass is less than 
the expansion rate during inflation. The moduli obtain a mass ~ (9(TeV) from 
supersymmetry breaking. They start oscillating with a large amplitude, possibly as 
big as Mp, when the Hubble parameter drops below their mass. Since moduli are 
only gravitationally coupled to other fields, their oscillations dominate the Universe 
while they decay very late. The resulting reheat temperature is below MeV, and is 
too low to yield a successful primordial nucleosynthesis. 

However, in our case iJjnf ^ TeV . This implies that quantum fluctuations 
cannot displace the moduli from their true minima during the inflationary epoch 
driven by MSSM flat directions. Moreover, any oscillations of the moduli will be 
exponentially damped during the inflationary epoch. Therefore our model is free 
from the infamous moduli problem . 



7. Cold dark matter and MSSM inflation 

Since is related to the scalar masses, sleptons {LLe direction) and squarks {udd 
direction), the bound on will be translated into the bounds on these scalar masses 



which are expressed in terms of the model parameters P5 |. 

The models of mSUGRA depend only on four parameters and one sign. These 
are uiq (the universal scalar soft breaking mass at the GUT scale Mq); mi/2 (the 



18t 



Note that the conditions in Eqs. (|6.6,6.7) can be simuhaneously satisfied easily. 
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universal gaugino soft breaking mass at Mq); Ao (the universal trilinear soft breaking 
mass at Mq) tan/5 = {H2){Hi) at the electroweak scale (where H2 gives rise to 
u quark masses and Hi to d quark and lepton masses); and the sign of /i, the 
Higgs mixing parameter in the superpotential (W^ = HH1H2). Unification of gauge 
couplings within supersymmetry suggests that Mq ^ 2 x 10^^ GeV. The model 
parameters are already significantly constrained by different experimental results. 
Most important constraints are: 



The hght Higgs mass bound of M^o > 114.0 GeV from LEP |64 . 



The 6 S7 branching ratio [|10|]: 2.2 x 10"^ < B{B X,7) < 4.5 x 10" 



In mSUGRA the x? is the candidate for CDM. The 2a bound from the WMAP 
PI gives a rehc density bound for CDM to be 0.095 < ^cdm/^^ < 0.129. 



The bound on the lightest chargino mass of M-± > 104 GeV from LEP 103 . 

Xi 



The possible 3.3 a deviation (using e+e" data to calculate the leading order 
hadronic contribution) from the SM expectation of the anomalous muon mag- 
netic moment from the muon g — 2 collaboration [p.04| . 



The allowed mSUGRA parameter space, at present, has mostly three distinct 
regions: (i) the stau-neutralino (fi — Xo), coannihilation region where Xo is the 
lightest SUSY particle (LSP), (ii) the Xo having a dominant Higgsino component 
(focus point) and (iii) the scalar Higgs (A°, H^) annihilation funnel (2M^i ~ M 



These three regions have been selected out by the CDM constraint. There stills exists 
a bulk region where none of these above properties is observed, but this region is 
now very small due to the existence of other experimental bounds. After considering 
all these bounds we will show that there exists an interesting overlap between the 



constraints from inflation and the CDM abundance [ ^ . 

We calculate at 0o and 0o is 10^^ GeV which is two orders of magnitude 
below the GUT scale. From this m^, we determine mg and mi/2 by solving the RGEs 



19 



The relationship between the two A terms, the trihnear, Aq and the non-renorniahzable A term 



in Eq.(9.17) can be related to each other, however, that depends on the SUSY breaking sector. For 



a Polonyi model, they are given by: A = (3 - \/3)/(6 - \/i)Aa Esf 
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Figure 5: The contours for different values of Ug and Sh are shown in the mo — 1^1/2 plane 
for tan f3 = 10. We used A = 1 for the contours. We show the dark matter allowed region 
narrow blue corridor, (g-2)^ region (light blue) for < 11 x 10~®, Higgs mass < 114 GeV 
(pink region) and LEPII bounds on SUSY masses (red). We also show the the dark matter 
detection rate by vertical blue lines. 

for fixed values of Aq and tan j3. The RGEs for are 

drr^ = Zl-{^Mlgl + -M^,) , (for LLe) 

dfji^ 1 o 

= ^(4M|^?I + , (forudd) . (7.1) 

Ml, M2 and M3 are U{1), SU{2) and SU{?>) gaugino masses respectively. 

After we determine mo and mi/2 from m^, we can determine the allowed values 
of from the experimental bounds on the mSUGRA parameters space. In order 
to obtain the constraint on the mSUGRA parameter space, we calculate the SUSY 
particle masses by solving the RGEs at the weak scale using four parameters of the 
mSUGRA model and then use these masses to calculate Higgs mass, BR\b 57], 
dark matter content etc. 
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Figure 6: The contours for different values of Ug and Sh are shown in the mo — 1^1/2 plane 
for tan f3 = 40. We used A = 1 for the contours. We show the dark matter allowed region 
narrow blue corridor, (g-2)^ region (light blue) for < 11 x 10"*^, 6 ^ 57 allowed region 
(brick) and LEPII bounds on SUSY masses (red). 

We show that the mSUGRA parameter space in figures ^ ^ for tan (3 = 10 and 

40 with the udd fiat direction using A = 1 In the figures, we show contours 

correspond to ris = 1 for the maximum value of 6h = 2.03 x 10~^ (at 2cr level) and 

Ug = 1.0, 0.98, 0.96 for 6h = 1-91 x 10~^. The constraints on the parameter space 

arising from the inflation appearing to be consistent with the constraints arising from 

the dark matter content of the universe and other experimental results. We find that 

tan/3 needs to be smaller to allow for smaller values of < 1. It is also interesting to 

note that the allowed region of m^, as required by the inflation data for A = 1 lies in 

the stau-neutralino coannihilation region which requires smaller values of the SUSY 

particle masses. The SUSY particles in this parameter space are, therefore, within 

the reach of the LHC very quickly. The detection of the region at the LHC has been 

^"^We have a similar figure for the flat direction LLe which we do not show in this paper. All the 
figures are for udd flat direction as an inflaton. 
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Figure 7: The contours for different values of Ug and 6h are sliown in the mo — ?7ii/2 plane 
for tan f3 = 10. We used A = 0.1 for the contours. We show the dark matter allowed region 
narrow blue corridor, g-2 region (light blue) for < 11 x 10~^, Higgs mass < 114 GeV 
(pink region) and LEPII bounds on SUSY masses (red). The black region is not allowed 
by radiative electroweak symmetry breaking. We use mt = 172.7 GeV for this graph. 



considered in refs ||105|| . From the figures, one can also find that as tan/? increases, 
the inflation data along with the dark matter, rare decay and Higgs mass constraint 
allow smaller ranges of mi/2- For example, the allowed ranges of gluino masses are 



765 GeV-2.1 TeV and 900 GeV-1.7 TeV for tan/? = 10 and 40 respectively |8 



So far we have chosen A = 1. Now if A is small e.g., A <, 10~^, we find that the 
allowed values of to be large. In this case the dark matter allowed region requires 
the lightest neutralino to have larger Higgsino component in the mSUGRA model. 
As we will see shortly, this small value of A is accommodated in 5*0(10) type model. 
In figure |^, we show Ug = 1, 0.98 contours for 6h = 1-91 x 10^^ in the mSUGRA 
parameter space for tan/3 = 10. In this figure, we find that rig can not smaller than 
0.97, but if we lower A which will demand larger and therefore Ug can be lowered 
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Figure 8: Contours of A for 6h = 1-91 x 10^^ in the ng-m^ plane. The blue band on the 
left is due to the stau-neutralino coannihilation region for tan (3 = W and the blue band 
on the right (which continues beyond the plotting range) denotes the focus point region. 

down to 0.92. 

In figure ^, we show the contours of A for different values of which are allowed 
by Hg and 6h = 1-91 x 10~^. The blue bands show the dark matter allowed regions 
for tan/? = 10. The band on the left is due to the stau-neutralino coannihilation 
region allowed by other constraints and the allowed values of A are 0.3-1. The first two 
generation squarks masses are 690 GeV and 1.9 TeV for the minimum and maximum 
values of allowed by the dark matter and other constraints. The gluino masses 
for these are 765 GeV and 2.1 TeV respectively. The band is slightly curved due to 
the shifting of 0o as a function A. (We solve for SUSY parameters from the infiaton 
mass at ^o). The band on the right which continues beyond the plotting range of 
the figure § is due to the Higgsino dominated dark matter. We find that A is mostly 
< 0.1 in this region and > 1.9 TeV. In this case the squark masses are much 
larger than the gluino mass since mo is much larger than mi/2- 
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8. Grand unified Models and Inclusion of Right-Handed Neu- 
trinos 



8.1 Embedding MSSM inflation in SU{5) or SO(IO) GUT 

As we have pointed out, mSUGRA makes a mild assumption that there exists a GUT 
physics which encompasses MSSM beyond the unification scale Mq Here we wish 
to understand how such embedding would affect infiationary scenario, for instance, 
would it be possible to single out either LLe or udd as a candidate for the MSSM 
inflaton. 

The lowest order non-renormalizable superpotential terms which lift LLe and 
udd are (see Eq. (|9.1D): 

{LLer {uddf 
Ml ' M3 • ' 

It is generically believed that gravity breaks global symmetries. Then all gauge 
invariant terms which are Mp suppressed should appear with A ~ Obviously 
the above terms in Eq. ( |8.1| ) are invariant under the SM. Once the SM is embedded 
within a GUT at the scale Mq, where gauge couplings are unified, the gauge group 
will be enlarged. Then the question arises whether such terms in Eq. ( |8.1| ) are 
invariant under the GUT gauge group or not. Note that a GUT singlet is also a 
singlet under the SM, however, the vice versa is not correct. To answer this question, 
let us consider SU{5) and 5*0(10) models separately. 

• SU{5y. 

We briefiy recollect representations of matter fields in this case: L and d belong 
to 5, while e and u belong to 10 of SU{5) group. Thus under SU{5) the 
superpotential terms in Eq. (OF) read 



5x5x10x5x5x10 
Ml • 

This product clearly includes a SU{5) singlet. Therefore in the case of SU{5), 
we expect that Mp suppressed terms as in Eq. (pTT|) appear with A ~ 0{1) 



^^We remind the readers that inflation occurs around a flat direction VEV (j>o ~ 10^^ GeV. Since 
00 ^ Mg, heavy GUT degrees of freedom play no role in the dynamics of MSSM inflation, and 
hence they can be ignored. 

^^If we were to obtain the [LLeY term by integrating out the heavy fields of the SU{5) GUT, 
then A = 0. This is due to the fact that SU{5) preserves B — L. 
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50(10): 

In this case all matter fields of one generation are included in the spinorial 
representation 16 of SO{10). Hence the superpotential terms in Eq. ( |8.1j ) 
are [16]^ under SO{10), which does not provide a singlet. A gauge invariant 
operator will be obtained by multiplying with a 126-plet Higgs. This implies 
that in 5*0(10) the lowest order gauge invariant superpotential term with 6 
matter fields arises at n = 7 level: 

16 X 16 X 16 X 16 X 16 X 16 X 126^ 



.3) 



Once 126h acquires a VEV, 50(10) can break down to a lower ranked subgroup, 
for instance SU (5). This will induce an effective n = 6 non-renormalizable term 
as in Eq. (|9J] ) with 

{126;^^0(Mo^ (84) 
Mp Mp ■ ^ ■ ^ 

Hence, in the case of 50(10), we can expect A ~ O{10~^ — 10^^) depending on 
the scale where SO (10) gets broken. 

We conclude that embedding MSSM in 50(10) naturally implies A ^ 1. Hence 
an experimental confirmation of the focus point region may be considered as an 
indication for 50(10). More precise determination of the spectral index Ug from 
future experiments (such as PLANCK) can in addition shed light on the scale of 
50(10) breaking. Smaller values of Ug (within the range 0.92 < < 1) point to 
smaller A, as can be seen from figure 6. This, according to Eq. ( ^.4|) , implies a scale 
of 50(10) breaking, i.e. (126//), which is closer to the GUT scale. 

Further note that embedding the MSSM within 50(10) also provides an advan- 
tage for obtaining a right handed neutrino. 

8.2 Including Right-Handed Majorana Neutrinos 

Eventually one would need to supplement MSSM with additional ingredients to ex- 
plain the tiny neutrino masses. Here we consider the most popular framework; the 
see-saw mechanism which invokes MSSM plus three RH (s)neutrinos A^i, A^2) A^s with 
respective Majorana masses Mj. By adding new superfields to MSSM, one can write 
a larger number of non-renormalizable gauge- invariant terms of the form in Eq. (|9.1| ). 
As a result, a given flat direction might be lifted at a a different superpotential level. 
Then a natural question arises that whether/how adding new superfields will affect 
the inflaton candidates, i.e. LLe and udd flat directions. 
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Since, Ni, 1 < i < 3, are SM singlets, we can write the following n = 4 superpo- 
tential terms: 

NiLLe Niudd , , 

Mp ' Mp ■ ^ ■ 

Note that these terms are also singlet under SU (5) and 50 (10). In the case of SU (5), 

the terms in Eq. (|8.5|) read 5x5x10x1, which includes a singlet. While in the case 

of S'O(IO), since belongs to the 16, the terms in Eq. ( p.5|) read 16 x 16 x 16 x 16, 

which includes a singlet. Hence both terms in Eq. (|8.5|) are allowed in SU{5) or 

50(10) embedding of MSSM as well ^3. 

We now analyze the case for two flat directions separately. 

• LLe: 

First let us consider the LLe fiat direction. Taking into account of the family 
indices, there are 5 independent D-fiat directions as such Within MSSM, 
there are three directions which are F-fiat at the n = 3 level, one of which 
survives until n = 6. However the term in Eq. (|8.5|) leads to three additional 
F-term constraints F^^ = 0, which are more than sufficient to lift the remaining 
direction at the n = 4 superpotential level 

Generically in this case we would expect LLe to be lifted by a non-renormalizable 
operator n < 6. 

• udd: 

Next consider the udd direction. With family indices taken into account, there 
are 9 independent D-fiat directions as such Within MSSM, 3 directions 



are lifted by n = 4 terms uude/Mp, while the remaining 6 will be lifted at the 
n = 6 level. Note that the superpotential term in Eq. (|8.5|) lead to three F- 
term constraints at the n = 4 level. Nevertheless, 3 directions will still survive 
until n = 6. 

Based on the above analysis, if we include the RH neutrinos, we conclude that 
udd direction is a more promising infiaton candidate than LLe. The reason is that 
the flatness of the former will not be lifted in the presence of physically motivated 
right handed neutrino fields in addition to that of the MSSM fields. 



^^In the case of 5*0(10) one can naturally obtain a right-handed neutrino. 

^^The gauge invariant LLe direction will survive until n = 6 if all Mi 3> ^o- However this is not 
a phenomenologically viable situation. 
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9. Few more examples 



Within MSSM there are other interesting possibihties of inflation which we will dis- 
cuss in this section. In the flrst section we will discuss inflation with Dirac neutrinos. 

9.1 Inflation with Dirac neutrinos 

As we know by now that the inflaton potential has to be cosmologically flat, which 
is suggestive of either a symmetry or a small coupling, or both. We will see that 
this property of the inflaton may be related to the smallness of neutrino masses. 
Identifying such a connection could have important ramiflcations and could lead to a 
more fundamental theory. The model we will use as an example will contain nothing 
but the MSSM and the right-handed neutrinos. We will show that a viable inflation 
in this model favors the correct scale for the neutrino masses. 

Let us consider the MSSM with three additional flelds, namely the right-handed 
(RH) neutrino supermultiplets. The relevant part of the superpotential is 



= W^MSSM + hNH„L. (9.1) 

Here N, L and are superflelds containing the RH neutrinos, left-handed (LH) 
leptons and the Higgs which gives mass to the up-type quarks, respectively. For con- 
ciseness we have omitted the generation indices. We note that the RH (s) neutrinos 
are singlets under the standard model (SM) gauge group. However in many exten- 
sions of the SM they can transform non-trivially under the action of a larger gauge 
group. The simplest example is extending the SM gauge group to SU (3)c x S'f/(2)w x 
U{1)y X [/(1)b-l, which is a subgroup of SO{10). Here B and L denote the baryon 
and lepton numbers, respectively. This is the model we consider here. In particular, 
the U{1)b-l prohibits the RH Majorana masses 

The active neutrino masses that arise from this are given by the usual seesaw 
relation {Hu)"^ / M | |106| , [L07| |, where {Hu) is the Higgs vacuum expectation value 
(VEV). Although the seesaw mechanism allows for small neutrino masses in the 
presence of large Yukawa couplings, it does not require the Yukawa couplings to be of 
order one. It still allows one to choose between the large Yukawa couplings and large 
Majorana masses on the one hand, and the small Yukawas, small Majorana masses. 



^^The monomials with B — L ~ will be also D-flat under U (1)b-l, while those with B — L ^ Q 
must be multiplied by an appropriate number of N superfields. In particular, NH„L is now a 
D-flat direction. 
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on the other hand. Viable models for neutrino mass matrices have been constructed 
in both limits, including the low-scale seesaw models [|108| , |109|| , in which the Yukawa 
couphngs are typically of the order of 

h ~ 10"^^ (9.2) 

or the same order of magnitude, as it would have in the case of Dirac neutrinos 
in order to explain the mass scale ~ (9(0.1 eV) corresponding to the atmospheric 
neutrino oscillations detected by Super-Kamiokande experiment. 

Let us now work in the basis where neutrino masses are diagonalized. There is 
a flat direction NsH^La spanned by the VEV of the lower and upper weak isospin 
components of H„ and L3, respectively. The scalar field corresponding to the fiat 
direction is denoted by 

,.«i±^. ,8.3) 

where the superscripts refer to the weak isospin components. One must now include 
the soft SUSY breaking terms, such as the mass terms and the A-teim. The A-terms 
are known to play an important role in Affleck-Dine baryogenesis [^, as well as in 
the inflation models based on supersymmetry [ |63|1 . 



The potential along the flat direction is found to be 



h'-^ Ah 

V{<j>) = -fct^' + + 0|cos {e + e, + Oa) <P' , (9.4) 

where the flat direction mass is given in terms of the soft masses of N^, Hu, and L3: 
m| = (^rn^fj^ + nij^^ + m| j /3. Here we have used the radial and angular components 
of the flat direction (pji + icpj = ^/2(f) exp (i6) , and 9h, Oa are the phases of the Yukawa 
coupling /13 and the A-term, respectively. We note that the above potential does not 
contain any non-renormalizable term at all. 

The last term on the right-hand side of eq. ( |9.4|) is minimized when cos {0 + 9h + 9a) 
— 1. Along this direction, V{(j)) has the global minimum at = and a local mini- 
mum at 00 ~ rn^/hs, as long as 

4m^ <A< 3V2m^ . (9.5) 



When the inequality in eq. (|9.5| ) is saturated, i.e., when A = Am^, then both flrst 
and second derivatives of V vanish at 0o, ^'(0o) = ^"(0o) = 0, and the potential 
becomes extremely flat in the radial direction, see Fig. H. We note that individually 
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9„ 



9 



Figure 9: The inflaton potential. The potential is flat near the saddle point where inflation 
occurs. 



none of the terms in eq. ( |9.4|) could have driven a successful inflation at VEVs lower 
than Mp = 2.4 x 10^^ GeV. However the combined effect of all the terms leads to a 
successful inflation without the graceful exit problem. 

Around 0o the field is stuck in a plateau with potential energy 



4hr 



(9.6) 



The first and second derivatives of the potential vanish, while the third deriva- 
tive does not. Around = 0o one can expand the potential as V{(f)) = V"(0o) + 
(l/3!)r"(0o)(0-0o)^ where 



V" 



(9.7) 



Hence, in the range [0o- A0, 0o + where A0 ~ Hl^/V'i^o) ~ {(pl/M^) > Hi^i, 
the potential is flat along the real direction of the inflaton. Inflation occurs along 
this flat direction. 

If the initial conditions are such that the flat direction starts in the vicinity of 
(j)Q with 0^0, then a sufficiently large number of e-foldings of inflation can be 
generated. Around the saddle point, due to the random fluctuations of the massless 
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field, the quantum diffusion is stronger than the classical force, Hi^ijlix > 0/i/inf 
for 

(00-0) / /m^02 1/2 _ / 3^3 xl/2 



At later times, the evolution is determined by the usual slow roll. The equation of 
motion for the field in the slow-roll approximation is 3H(p = — (1/2) V"(0o)(0 — 



A rough estimate of the number of e-foldings is then given by 



Hd(j) / 



where we have assumed V(0) ~ (0 — 0o)^V""'(0o) (this is justified since V^'(0o) and 
V"{(f)o) are both small). We note that the initial displacement from 0o cannot be 
much smaller than -ffjnf, due to the uncertainty from quantum fiuctuations. 

Infiation ends when the slow roll parameters become ~ 1. It turns out that 
|?7| ~ 1 gives the dominant condition 

.3 

(9.10) 



The total number of e-foldings can be computed as [0, ^ 



^e-(A^y^=(^] , (9.11) 



\m^Mp J \hlMp 
evaluated after the end of diffusion, see eq. (|9.8| ), when the slow-roll regime is 
achieved. 



Let us now consider adiabatic density perturbations. As in Ref . ^ ^ , one 
finds 

1 Hlj hlMp ^ 

^ 3 ~ A/^BE • (9-12) 

Stt 3m^ 

In the above expression we have used the slow roll approximation ^ —V"'{(f)Q){(f)Q — 
0)2/3ifinf, and eq. (|1Q.16|) . The exact number depends on the scale of infiation and 
on when the Universe becomes radiation dominated (we note that full thermalization 
is not necessary as it is the relativistic equation of state which matters). In our case 
A/coBE < 60 as we shall see below. 

The spectral tilt of the power spectrum and its running are 

4 

Us = l + 2r]-Qe^l- — , (9.13) 

A/cOBE 

^ (9.14) 



dink AT^obe' 
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cf. 0. (We note that e < 1 while r] = -2/A/'cobe.) 

It is a remarkable feature of the model that for the weak-scale supersymmetry 
and for the correct value of the Yukawa coupling, namely, 



~ 100 GeV - 10 TeV , ~ 10"^^ , (9.15) 

the flat direction NsH^La leads to a successful low scale inflation near 0o ~ (lO^'^ — 10^^) GeV <C 
Mp, with 

V ~ 10^2 _ 1036 Gev^ , ffinf ~ 10 MeV - IGeV , 

K ~ 10=^ , Tma. ~ 10' - 10' GeV . (9.16) 

The total number of e- foldings driven by the slow roll inflation, A/'e ~ 10^, is more 
than sufficient to produce a patch of the Universe with no dangerous relics. Those 
domains that are initially closer to 0o enter self-reproduction in eternal inflation. 
Since the infiaton, NaH^Ls, couples directly to MSSM particles, after inflation the 
field oscillates and decays to relativistic MSSM degrees of freedom. The highest 
temperature during reheating is T^ax ~ V^^^. This temperature determines the 
total number of e-foldings required for the relevant perturbations to leave the Hubble 
radius during infiation; in our case it is roughly A/cobe ~ 50. 

Despite the low scale of infiation, the fiat direction can generate density perturba- 
tions of the correct size for the parameters listed above. Indeed, from eqs. ( |9.26|J9.29| ), 
and (|9.15|), we obtain: 6h ~ 10"^. Following the discussion of Section 4.4, we obtain 



0.91 < Us < 1.0 with a negligible running. The spectral tilt agrees with the current 



WMAP 3- years' data within 2a ||T0[ . The tensor modes are negligible because of the 
low scale of infiation. 

We emphasize that the VEV of the fiat direction is related to the Yukawa coupling 
that can generate the Dirac neutrino mass ~ 0.1 eV. The scale of the neutrino mass 
appears to be just right to get the correct amplitude in the CMB perturbations. 

The infiaton has gauge couplings to the electroweak and U{1)b-l gauge/gaugino 
fields. It therefore induces a VEV-dependent mass ~ g{(f)) for these fields {g denotes 
a typical gauge coupling). After the end of infiation, starts oscillating around the 
global minimum at the origin with a frequency m,^ ~ lO'^iJinf, see eq. (|9.16| ). When 
the infiaton passes through the minimum, (0) = 0, the induced mass undergoes non- 
adiabatic time variation. This results in non-perturbative particle production |^ . 
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9.2 Inflation in gauge mediated scenarios 

In a Gauge Mediated Supersymmetry Breaking (GMSB) scenario the two-loop correc- 
tion to tlie fiat direction potential results in a logarithmic term above the messenger 
scale, i.e. > Ms ||110|| . Together with the A-term this leads to the scalar potential 



where Mp ~ {msusY x ^sY^"^ and msusY ~ 1 TeV is the soft SUSY breaking mass 
at the weak scale. For cf) > M|./m3/2, usually the gravity mediated contribution, 
7713^20^) dominates the potential where 1713/2 is the gravitino mass. Here we will 
concentrate on the VEVs Mg <^ (f> <, Mp / 1713/2. 

Although individual terms are unable to support a sub-Planckian VEV inflation. 



but as shown in Refs. Q , a successful inflation can be obtained near the 

saddle point, which we flnd by solving, V'{(l)o) = V"{(j)o) = (where derivative is 
w.r.t 6). 



n-3 /i^2 I \ V("-l) 

(9.18) 



A„ y{n-l){n-2) 
A=^±^^ro-'. (9.19) 

In the vicinity of the saddle point, we obtain the total energy density and the third 
derivative of the potential to be |^ : 



^(4) = M^p 



In 



bl\ 3n-2 



MlJ n{n 



(9.20) 



V"\<Po) = 4n(n - l)M^0o-^ . (9.21) 

There are couple of interesting points, flrst of all note that the scale of inflation is 
extremely low in our case, barring some small coefficients of order one, the Hubble 
scale during inflation is given by: 

i/inf ~ M|/Mp ~ 10"^ - 10"^ eV , (9.22) 

for Mp ~ 1 — 10 TeV. For such a low scale inflation usually it is extremely hard to 
obtain the right phenomenology. But there are obvious advantages of having a low 
scale inflation, Mp ^ -f^inf- The supergravity corrections and the Trans-Planckian 
corrections are all negligible [KH], therefore the model predictions are trustworthy. 
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Perturbations which are relevant for the COBE normahzation are generated a 
number A/cobe e-foldings before the end of inflation. The value of A/cobe depends 
on thermal history of the universe and the total energy density stored in the infia- 
ton, which in our case is bounded by, Vq < 10^^ (GeV)^. The required number of 
e-foldings yields in our case, A/cobe ~ 40 [^], provided the universe thermahzes 
within one Hubble time. Although within SUSY thermalization time scale is typi- 
cally very long [|^, however, in this particular case it is possible to obtain a rapid 
thermalization. 

Near the vicinity of the saddle point, ^o, the potential is extremely fiat and one 
enters a regime of self-reproduction [Q. The self-reproduction regime lasts as long 
as the quantum diffusion is stronger than the classical drag; Hinf/2TT > (p/Hinf, for 
0s < < 00, where 0o ~ 0s — Mp {(po/ Mp)^^"^ . From then on, the evolution is 



governed by the classical slow roll. Inflation ends when e 
~ 5 where [^| 



1, which happens at 



(9.23) 



V2n{n - l)M|Mp ' 

Assuming that the classical motion is due to the third derivative of the potential. 



V'{4>) ^ (l/2)\/'"(0o)(0- 
period is found to be ^Tj : 



the total number of e-foldings during the slow roll 



This simplifies to E7 



Hd(p 



2Vo4>l 



An{n - 1)M|M| V0o - 



Kot ^ 0f /(M^/'Mp) 



(9.24) 



(9.25) 



Let us now consider the adiabatic density perturbations. Despite Hi^f <^ 1 eV, the 
fiat direction can generate adequate density perturbations as required to explain to 
match the observations. Recall that inflation is driven by V" 7^ 0, we obtain |^ 



(9.26) 



Note that for Mp 
10[| , when 



tions 



10 TeV, and Acobe ~ 40, we match the current observa- 
()o ~ 10^^ GeV. The validity of Eq. (pTzl) for such a large VEV 
requires that M| > (10^^ GeV) x m3/2. For Mp ~ 10 TeV this yields the bound 
on the gravitino mass, m^i2 < 1 MeV, which is compatible with the dark; matter 
constraints as we will see. 
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We can naturally satisfy Eq. ( |9.26| ) provided, n = 6. The non-renormalizable 
operator, n = 6, points towards two MSSM flat directions out of many, 

LLe and udd . (9.27) 



As we discussed before in |Q, these are the only directions which are suitable for 
inflation as they give rise to a non-vanishing A-term. Note that the inflatons are 
now the gauge invariant objects. The total number of e-foldings, during the slow roll 
inflation, after using Eq. ( |9.25D yields, 

Mtot ~ 10' . (9.28) 

While the spectral tilt and the running of the power spectrum are determined by 
A/boBE ~ 40 < Alot. 

4 

= l + 2r/-6e^l-— ^ 0.90, (9.29) 

AcOBE 

- 16er/ - 24e2 - 2^2 ~ --.^ 10"' 



(9.30) 



where = MpV'V" /V"^ . Note that the spectral tilt is slightly away from the 2a 
result of the current WMAP 3 years data, on the other hand running of the spectrum 



is well inside the current bounds [10 



At first instance one would discard the model just from the slight mismatch 
in the spectral tilt from the current observations. However note that our analysis 
strictly assumes that the slow roll infiation is driven by V"'{(f)Q). This is particularly 
correct if V^'(0o) = and V'(0o) = 0. Let us then study the case when V^'(0o) 7^ 0, 



as discussed in Ha, 50 



The latter case can be studied by parameterizing a small deviation from the 
exact saddle point condition by solving near the point of infiection, where we wish 
to solve V'^cpo) = and we get up-to 1st order in the deviation, 6 < 1, 

A = A{1- 6), 00 = 00 f 1 - -^^j) > (9.31) 

\ n{n — 2) / 

with A and 0o are the saddle point solutions. Then the 1st derivative is given by 

Ti 1 

V^'(0o) = 4 (9.32) 

n — 2 

Therefore the slope of the potential is determined by, V'{(f)) ~ V'{(f)o) + {l/2)V"'{(f)o){(f)— 
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Note that both the terms on the right-hand side are positive. The fact that 
^'(0o) 7^ can lead to an interesting changes from the saddle point behavior, for 
instance the total number of e-foldings is now given by 



tot 



Ml 



(f>Q 



(9.33) 



First of all note that by including V, we are slightly away from the saddle point 
and rather close to the point of inflection. This affects the total number of e-foldings 
during the slow roll. It is now much less than that of A/'tot, i-e. A/'tot ^ 10^, see 
Eq. dpI). 

When both the terms in the denominator of the integrand contributes equally 
then there exists an interesting window. 



k/8<6 < k/2. 



where 



n-2 



K = 



n[n 



1)^ 



In 



Ml 



3n-2 



n[n 



MM 



(9.34) 



(9.35) 



'P-'^COBE 

The lower limit in Eq. ( |9.34|) is saturated when V(0o) = 0, while the upper limit is 
saturated when A/'tot — Acobe — 40. It is also easy to check that there will be no 
self-reproduction regime for the field values determined by 5. 

It is a straightforward but a tedious exercise to demonstrate that when the upper 
limit of Eq. ( |9.34 ) is saturated the spectral tilt becomes ^ 1, when the lower limit 
is satisfied we recover the previous result with Us = 0.9. This value, Ug —>■ 1, can be 
easily understood as 0cobe 0o (where 0cobe corresponds to the VeV where the 
end of inflation corresponds to Acobe ~ 40), in which case, r/ — > 0. Therefore the 
spectral tilt becomes nearly scale invariant. We therefore find a range [HR BUI E7|, 



0.90 <n,< 1 



(9.36) 



whose width is within the 2a error of the central limit |jTO[ . Similarly the running of 
the spectral tilt gets modified too but remains within the observable limit while 



the amplitude of the power spectrum is least affected 45, 50, 47, 48 



Let us now discuss the issue of reheating and thermalization. Important point 
is to realize that the infiaton belongs to the MSSM, i.e. LLe and udd, both carry 



^^A similar exercise can be done for the running of the spectral tilt and the running lies between 
-le/A/^oBE < dnjdlnk < -A/M^^be ® 0- 
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MSSM charges and both have gauge couphngs to gauge bosons and gauginos. After 
inflation the condensate starts oscillating. The effective frequency of the inflaton 
oscillations in the Logarithmic potential, Eq. ( p.l7| ), is of the order of Mp/cpo, while 



the expansion rate is given by i^inf ~ Mp/Mp. This means that within one Hubble 
time the inflaton oscillates nearly Mp/0o ~ 10"^ times. The motion of the inflaton is 
strictly one dimensional from the very beginning. During inflation, the imaginary 
direction is very heavy and settles down in the minimum of the potential. This also 
prohibits fragmentation of the flat direction to form Q-balls ||TP|, Although one 



could still argue positive and negative charged Q-ball formation ||101|| , but they do 
not form adequately to alter reheating and thermalization in any signiflcant way. 

An efficient bout of particle creation occurs when the inflaton crosses the origin, 
which happens twice in every oscillation. The reason is that the fields which are 
coupled to the infiaton are massless near the point of enhanced symmetry. Mainly 
electroweak gauge fields and gauginos are then created as they have the largest 
coupling to the fiat direction. The production takes place in a short interval. Once 
the infiaton has passed by the origin, the gauge bosons/gauginos become heavy by 
virtue of VeV dependent masses and they eventually decay into particles sparticles, 
which creates the relativistic thermal bath. This is so-called instant preheating 



mechanism ||9^. In a favorable condition, the fiat direction VeV coupled very weakly 
to the fiat direction infiaton could also enhance the perturbative decay rate of the 
infiaton |]^. In any case there is no non-thermal gravitino production as the 



energy density stored in the infiaton oscillations is too low. 

A full thermal equilibrium is reached when a) kinetic and b) chemical equilibrium 
are established The maximum temperature of the plasma is given by 

Tr ~ [V{<f)o)]'^^ --Mf<, lOTeV , (9.37) 

when the fiat direction, either LLe or udd evaporates completely. This naturally 
happens at the weak scale. There are two very important consequences which we 
summarize below. 

9.2.1 Cold electroweak Baryogenesis 

The model strongly points towards electroweak baryogenesis within MSSM. Note 
that the reheat temperature is sufficient enough for a thermal electroweak baryoge- 
nesis IP-Ill . 
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However, if the thermal electroweak baryogenesis is not triggered, then cold elec- 
troweak baryogenesis is still an option ||112|| . During the cold electroweak baryogene- 



sis, the large gauge field fiuctuations give rise to a non-thermal sphaleron transition. 
In our case it is possible to excite the gauge fields of SU{2)l x U{1)y during instant 
preheating provided the infiaton is LLe. The LLe as an infiaton carries the same 
quantum number which has a, B — L anomaly and large gauge field excitations can 
lead to non-thermal sphaleron transition to facilitate baryogenesis within MSSM. 

9.2.2 Gravitino dark matter 

Within GMSB gravitinos are the LSP and if the i?-parity is conserved then they are 
an excellent candidate for the dark matter. There are various sources of gravitino 
production in the early universe ||113| , |114|| . However in our case the thermal produc- 



tion is the dominant one and mainly helicity ±1/2 gravitinos are created. Gravitinos 
thus produced have the correct dark matter abundance for |ll(J, 



r^V^V (9 38) 

100 KeV few VI TeV/ VI TeV/ ' ^' ' 



where Mg is the gluino mass. For ms/2 ^ 100 KeV, Eq. ( |9.38|) is easily satisfied for 
Mg ~ 1 TeV and Tr ^ 10 TeV. We remind that for O(keV) ^ mg/s < 100 KeV 



gravitinos produced from the sfermion decays overdose the universe [IIC . 



10. Quantum initial conditions, cosmological constant and a 
string landscape 

One important result is that the MSSM infiation produces many e-foldings of slow 
roll infiation, AC ~ 10^, with a preceding self-reproduction regime. One difficulty of 
MSSM infiation, though, is that it requires a fine-tuning of the initial conditions for 



slow-roll |H5 



Our aim is to address this initial condition issue by studying the quantum fiuc- 
tuations of the MSSM infiaton during a false vacuum infiation which naturally fine 
tunes the initial conditions for MSSM infiation, i.e. why the MSSM infiation occurs 
near the saddle point of the potential. The idea is to have an early bout of false 
vacuum during which the MSSM infiaton will obtain false vacuum induced quantum 
fiuctuations from the point of enhanced symmetry to the saddle point. Once the 
false vacuum tunnels to a lower cosmological constant or the vacuum energy density 
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comparable to that of the MSSM inflation then the latter potential takes over for a 
last stage of MSSM driven inflation. 

Without a prior phase of inflation it would be difficult to explain the initial 
condition otherwise. For example, one may not invoke thermal effects to trap the 
flat direction near the saddle point as it does not correspond to a false minimum, or 
a point of enhanced symmetry. As we shall see, old inflation and MSSM inflation 
complement each other nicely. 

Also, there is some evidence that string theory has a "landscape" of metastable 
vacua with varying cosmological constant, moduli VEVs (and therefore couplings), 
supersymmetry (SUSY) breaking scale, and so on, which can be studied using sta- 
tistical arguments (see ||115| , |116| , |23| for reviews). Since many of these metastable 



vacua have large energy densities (as we describe below), it seems likely that the 
early universe could have existed as a de Sitter spacetime with a large cosmological 
constant, which then decayed by tunnelling, as in old inflation [Jl|. In fact, such a 
cosmology with multiple stages of inflation ||75[ provides a mechanism by which the 
full landscape of vacua is populated, as in | |117| , |118| , |119| |. With some caveats, this 



mechanism can also relax the cosmological constant quickly ||117| , |118| , |119| , |120| , |121 
In particular, the cosmological constant must be able to decay even though other 
sectors will dominate once the energy density reaches ~ 1 (TeV)^. 

One obvious worry about this picture is that, if the universe were to tunnel 
out of the false vacuum, then the universe would be devoid of any entropy as the 
nucleated bubble would keep expanding forever with a negative spatial curvature. 
Such a universe would have no place in a real world, so it is important that the last 
stage of inflation be driven in an observable sector, such as in the case of MSSM 
inflation which not only produces the right amplitude of density spectrum but also 
provides us with a desired thermal entropy,and cold dark matter abundance. 

10.1 Old Inflation on the Landscape 

The most basic fact about the landscape (which is usually not discussed, being of little 
interest for present-day physics) is that the cosmological constant is generically large. 



A simple way to see that follows |p.l8|| , which describes the landscape contribution 
to the cosmological constant as arising from string theory flux. In this picture, the 
vacuum energy 



1^ = M|A = M|,Ao + ^ , (10.1) 
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where Cj < 1 are constants and rij are flux quantum numbers (note that A has 
dimension mass^ in our notation as it enters the Einstein equation as h^Qf^y). It is 
clear that large A corresponds to a large radius shell in the space of flux quanta, so 
larger A will have more possible states. All in all, string theory (from the landscape 
point of view) could have from 10^''° to even 10^°°° vacua ||1 1 5| , |116| , p3| , |118|| , with 



the vast majority of those having large cosmological constants. 
10.2 Decay time scales and inflation 

The decay rate per volume (by tunnelling) of a metastable vacuum to the nearest 
neighboring vacuum^'' takes the form 

r/1/ = Cexp(-A5£) , (10.2) 

where C is a one-loop determinant and ASe is the difference in Euclidean actions 
between the instanton and the background with larger cosmological constant. The 
determinant C can at most be C < M|,, simply because Mp is the largest scale 
available, and estimates (ignoring metric fluctuations) give a value as small as C ~ 



r~ , with r the instanton bubble radius |[L24|| . If we therefore look at a decay rate in 
a (comoving) Hubble volume, we find 

/\/f4 

r<-^exp(-A^^) . (10.3) 

Especially with a large Hubble scale, the associated decay time is much longer than 
1/H, given that typically /S.Se ^ 1- 

In fact, as given in ||125[| following [ |12(j| , |127| , |128| , the Euclidean action takes the 



form 

ASe = 27r , (10.4) 

where the bubble radius and effective tension Tg are given in terms of dimensionless 
cosmological constants A = MpA/r^, and the actual tension of the bubble wall r. 
The full formulae are listed in the appendix. 

Unfortunately, a given tunnelhng process in the landscape probably has A+ ~ 1 
and A_ < 1. However, it may be instructive to look at three limits, A+ 0, oo. In 



27 



To avoid subtleties, we require that the both states have nonnegative energy density. See, 



for example, |122| for issues in the negative A case. Also, [123 has discussed the importance of 
negative A vacua in possibly separating parts of the landscape from each other. We adopt the view, 
as discussed in that paper, that the landscape is sufficiently complicated that there are no isolated 
regions. 
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the latter case, we simultaneously take A_ — > or A_ ^ A+ (these are the same for 
A+ — > 0). As A+ vanishes, we find 

A^^ ^ 247c'^ = 247r2^ . (10.5) 
T A+ A+ 

(The tension surprisingly drops out!) For A+ —>■ cxd, we find the limits 

ASe ^ I2n^^ = I2n^^ (A_ ^ 0) , (10.6) 
r^A+ A+ 

ASe^QVStt^ 1^ = 673772—- (A_^A+). (10.7) 

^2^3/2 ^3/2 V 

The only parametrically different behavior is the last limit as A± ^ oo, which means 
either that the tension is very small or that the cosmological constants are both very 
large (or that Mp — > oo, corresponding to field theory without gravity). Technically, 
the instanton approximation can break down in the limit ( p.0.7| ), if the Euclidean 
action becomes much smaller than one (as a full quantum mechanical treatment 
would be necessary), but the important point is that F/if > 1 if that is the case. 
From ( |10.3|) , we might expect that T/H ^ 1, but we should remember that the 



coefficient of the exponential in ( |10.3| ) is an upper limit and that quantum mechanical 
effects could limit the decay rate to F ~ if. In section |10.3| , we give an improved 
form of the argument of ||1 18|| that shows how old inflation on the landscape can 
source MSSM inflation without inflaton trapping and with only small jumps in the 
cosmological constant. 

Note that, when MSSM inflation starts, the "bare" cosmological constant (that 
not associated with the MSSM inflaton) might still be considerably larger than the 
present value. This means that further instanton decays should take place to reduce 
the bare cosmological constant, and these decays should occur during MSSM inflation 
in order for the bubble regions to grow long enough. Since it seems likely that the 
instanton bubble tension will be large compared to the scale of MSSM inflation, the 
decay rate will be given by (|10.5|) , which is highly suppressed. This would then 
require MSSM inflation to last for extremely many e-foldings. Fortunately, MSSM 
inflation naturally includes a self-reproduction (eternal inflation) regime prior to 
slow-roll [0, In other words, it seems that a low-scale eternal inflation is a 

necessary part of using decays to solve the cosmological constant problem, and MSSM 
inflation includes it naturally. Eternal inflation solves several of problems relating 
false vacuum inflation to observations, and it is particularly viable at the TeV scale 
of MSSM inflation. 
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10.3 Initial quantum kicks for an MSSM fields 

Let us now discuss the what happens to the observable sector in the background of 
inflation spacetime. As we have noticed, the universe cascades from large cosmologi- 
cal constant to another somewhat smaller as time progresses. Given the decay time, 
it is fairly evident that huge number of e-foldings can be generated. 

During the false vacuum inflation the energy density of the universe, and hence 
the expansion rate -fffaise, remains constant for a given A. 

The flat direction potential receives corrections from soft SUSY breaking mass 
term, the non-renormalizable superpotential correction, and corrections due to the 
large Hubble expansion for a minimal choice of Kahler potential. The Hubble cor- 
rection is relevant only when <C H{t), which is parameterized by two constants 



c ~ and a ~ C(l), result in |6|, |6|, |T9|, 
1 



V=-{ml + cHl,J\<p\' + 



{A + aiffaisc)A^ ^^^_3 + h.c.; 



|0|2(n-l) 
Ml 



+ A^ '"' , (10.8) 



where the soft SUSY breaking mass term is generically small compared to the Hubble 
expansion rate of the false vacuum, ~ 1 TeV ^ i^faise- We deflne $ = 06*""^, and 
the dynamics of in an inflationary background depends on c and a. Therefore we 
consider different cases separately. 

It turns out that the dynamics of the MSSM flat direction largely follow the 



physics discussed in ||118|| . We review and elaborate on their discussion in the context 
of the MSSM and give an important improvement on their result in the context of 
negligible Hubble corrections. 

10.3.1 Positive Hubble induced corrections 

A positive Hubble induced correction provides c ~ +(9(1), a ~ 0{1). This is a 
typical scenario when the Kahler potential for the string modulus comes with a 
canonical kinetic term. Although this could be treated as a special point on the 
Kahler manifold, it is nevertheless important to discuss this situation. This is also 



•^*The origin of the Hubble induced terms is due to couplings between the modulus which drives 
the false vacuum inflation and the MSSM sector. Apriori, even in string theory, the Kahler potential 
for the modulus is not well known at all points of the parameter space. Similarly, within the MSSM, 
the Kahler potential is unknown. Therefore, the coefficients c, a are not fixed. For a no-scale type 
model, the Hubble induced corrections are vanishing at tree level; however, they do appear at one 
loop |c| - 10-2 g ||. 
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the simplest scenario out of all possibilities. A generic flat direction gets a large 
mass of the order of Hubble expansion rate; therefore, its fluctuations are unable to 
displace the flat direction from its global minimum. 

On phenomenological grounds, this is an uninteresting and undesirable case, 
since the bubble does not excite any of the MSSM fields. Therefore, the bubble 
remains empty and devoid of energy with no graceful exit of inflation from the 
false vacuum. The universe continues cascading to smaller A with smaller -fffaise- 
Eventually, when iffaise < 'fn^, the MSSM flat directions would be free to move. 
However, through this time the fields were never displaced from their minimum and 
therefore the dynamics of the flat directions would remain frozen. The universe would 
be cold, as it was before, and the spatial curvature would remain negative. It is fair 
to say that, on phenomenological grounds, such a universe is already ruled out. This 
paves the way for more interesting scenarios, which we discuss next. 

10.3.2 Negligible Hubble induced corrections 

In this case, the potential is not affected by the false vacuum inflation at all, namely 
|c|, \a\ <^ 1. So long as V"{(f)) <^ -f^Mse^ ^^e flat direction field makes a quantum 
jump of length iffaise/27r within each Hubble time.^^ These jumps superpose in 
random walk fashion resulting in 

cm) = (10.9) 

\ at / fluctuations 47r^ 

On the other hand, the classical slow roll due to the potential leads to 

('My) = (10.10) 

V (It ' slow roll "J-TZfalse 

For a massive scalar field V{(f)) ~ m|(/)^/2, the combined effects yield B 



2ml 
exp ( -— — t 

J -"false 



The maximum field value 



3 -^filse 



10.11) 



;io.i2) 



at which the slow roll motion ( |10.10|) counterbalances the random walk motion ( |10.9| ) 
is reached for At ^ 3iyfaise/2m^. This amounts to a number 

3 / -f^false^ 



Malsc»^P^ . (10.13) 



^^To be more precise, the quantum fluctuations of have a Gaussian distribution, and the r.m.s. of 
modes which exit the horizon within one Hubble time is i/faisc/27r. 
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of e-foldings of inflation in the false vacuum. 

In the absence of Hubble induced corrections, the potential in ( |10.8| ) has a saddle 
point at 



where V'{<Pq) = V'icpo) = while V"'{(po) + 0, provided that = 8{n - l)m^.30 
For ~ 100 GeV - 10 TeV and A„ ~ C(l), and for n = 6, the VEV is 0o ~ 
10^"^ — 10^^ GeV. The suitable flat directions are LLe and udd, which are lifted by 
n = 6 superpotential terms and also have a non-zero A-term as required by the 
condition for a saddle point ^.^^ 



For (f) < (pQ the mass term dominates the flat direction potential. Then quantum 
fluctuations can push to the vicinity of 0o if 4'r.m,.s > 0o-^^ This, according to 
(|10.12D, requires that 



i^fai.e > (^j (m^0o)'/' ^ 10' GeV. (10.15) 



The number of e-foldings needed for this to happen is 



2 



iV,.,«<^^j 10'^. (10.16) 

Indeed for iffaisc > 0o the inherent uncertainty due to quantum fluctuations implies 
that > 00 within one Hubble time. 

A last stage of MSSM inflation with an expansion rate 



■^"^This can happen in the gravity mediated case, where A ^ ~ W3/2, where ~ C(l TeV) 
is the gravitino mass. The situation is quite different in the case of a gauge mediated SUSY breaking 
scenario Q. 

■^^In order to have successful inflation, we require the condition — 8(n — l)m^ to be satisfied to 
one part in 10^. Although this requires a fine tuning, SUSY can allow it to be maintained order by 
order if Ajm^ acts as an infrared fixed point of the renornialization group flow [ p3[ . Also, as noted 
in Ref. this tuning can be explained naturally by the landscape picture. For larger deviations 



there is a point of inflection with large V{<^q) (or a negligible A-term discussed in |69|), or a pocket 
of false minimum Neither case leads to a slow roll inflation within the MSSM. 

■^^Due to the Gaussian distribution of fluctuations, the probability of having ^ (pr.m.s is expo- 
nentially suppressed. 
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starts if V{(f)) dominates the energy density of the universe, i.e. V{(t>) > SHf^^^^Mp. 
An observationally consistent inflation in the slow roll regime requires that the dis- 
placement from the saddle point satisfy |0 — 0o| < where H5 



6^ 

A0 = ^ , ^ 10^ GeV. (10.18) 

An{n - l)Ml ^ ' 

In the landscape, the universe can begin in a false vacuum with arbitrarily large -fffaise 
(as long as -f^faise ^ ^p)- Therefore, we generically expect that is quickly pushed 
to field values ^ 0o- However, iffaise slowly decreases as a result of tunnelling 
to vacua with smaller cosmological constant For a massive scalar field with the 
potential V^(0) ~ m|0^/2, this implies a gradual decrease of (pr.m.s, see ( |10.12|) . Since 
quantum fluctuations can at most push to (pr.m.s, this also implies that is slowly 
decreasing in time. Indeed for iffaise < 10^ GeV, we find that < 0o, irrespective 
of how large initially was. This is the case in the discussion of | |118| ], which is 
why ||118|| requires a jump from large cosmological constant directly to the slow-roll 
inflationary stage. 

Note, however, that the potential becomes very flat around 0o as a result of the 
interplay among different terms in ( |10.8| ). In fact, for |0 — 0o| ^ 0o we have p^ , ^ 



1 TTi^ 

V{<P) ^ y(0o) + ^-^(0 - 0o)' . (10.19) 

3{n - 2)^ 00 

Once -fffaise ~ 10^ GeV, reaches this plateau (from above). Within the plateau, 
V'{(f)) becomes increasingly negligible, and so does the classical slow roll; they exactly 
vanish at = 0o. In consequence, quantum jumps dominate the dynamics and freely 
move throughout the plateau. Hence the flatness of potential, which is required 
for a successful MSSM inflation, guarantees that will remain within the plateau 
during the landscape evolutionary phase. It is possible to generalize this argument 
to other models of low-scale inflation; we see that a sufficiently flat inflaton potential 
can trap inflaton fluctuations in the slow-roll region. 

The flat direction eventually dominates the energy density of the universe when 
-f^faise < 1 GeV, see ( |10.171 ). MSSM inflation then starts in those parts of the universe 



■^^Note that we need to stay in an MSSM-like vacuum all the way until MSSM inflation begins. 
Given the scarcity of MSSM-like vacua in the landscape, the probability of tunnelling from one 
such vacuum to another is ~ 10~^. However, due to eternal inflation in the false vacuum, the 
physical volume of the universe increases by a factor of exp(3iJfaisc/r) within a typical time scale 
for bubble nucleation (F is the false vacuum decay rate, see subsection 2.2). This easily wins over 
the suppression factor 10^^ for F < 9i?faisc- 
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which obey ( 10.18 ). Having </> so close to 0o is possible since the uncertainty due to 
quantum fluctuations is ~ 1 GeV at this time. 

The false vacuum inflation paves the way for an MSSM inflation inside the nu- 
cleated bubble. The modulus which was responsible for the false vacuum inflation 
continues tunnelling to minima with smaller (eventually the currently observed) cos- 
mological constant inside the MSSM inflating bubble. The modulus will oscillate 
around the minimum of its potential as the curvature of its potential dominates over 
the Hubble expansion rate. The fate of oscillations would depend on the coupling 
of the modulus to the MSSM flelds. For a coupling of gravitational strength, os- 
cillations are long-lived. Moreover, the decay will be kinematically forbidden if the 
decay products are coupled to the flat direction (which has obtained a large VEV 
^ 0o)- However, such details are largely irrelevant as the flat direction dominates 
the energy density of the universe and drives inflation, diluting the energy density in 
oscillations. 



In particular, we note that MSSM inflation has a self-reproducing regime ^ 
because V'{(f)) is extremely small and the potential becomes very flat close to (pQ. 
The observable part of our universe can spend an arbitrarily long time in the self- 
reproducting regime before moving into the standard slow-roll inflation. During this 
period, the cosmological constant can continue to decay and eventually settle at an 
observationally acceptable value. 

Before we conclude this discussion, let us remind the readers that obtaining 
the flat direction VEV, 0o, depends on the false vacuum inflation, which requires 
-f^faise > 10^ GeV at some time (see ( |10.15| )). This condition, although very probable 
in the landscape picture, need not be satisfied always. Then the quantum fluctuations 
would not be large enough to push the flat direction to the vicinity of 0o as required 
for a flnal stage of MSSM inflation. This would therefore ruin the inflationary and 
phenomenological predictions. This can be avoided if the flat direction is trapped in 
a false minimum which evolves with time, as discussed below. 

10.3.3 Negative Hubble induced corrections 

The case with c ~ —0{1) may arise naturally for non-minimal Kahler potential ^8 



65[| . For -fffaisc ^ the potential in ( p.0.8|) becomes tachyonic, and its true minimum 

\n\/ In — 2 / 



is located at 

l/(n-2) 

|C 
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Note that 0min is initially larger than ^o, see (|10.14|) . 

The curvature of the potential at the minimum is V^"(0min) = (4^~7)|c|iJ£^jgp ^ 
-^faise- This implies that is driven away from the origin, due to quantum fluctua- 
tions, and quickly settles down at 0min- The MSSM fiat direction is trapped inside the 
minimum, held due to false vacuum inflation, and gradually tracks the instantaneous 
value of 0inin as -fffaise dccreases.^"^ 

The minima at = and = 0min become degenerate when iffaise ~ "^0- For 
-f^faise the true minimum is at = and the one at 0min will be false. The 



Hubble induced corrections are subdominant in this case. We then find 



\/"(0^i,) ~ 2{n - 2) v^/ffaisem,^ . (10.22) 



The field can track down the instantaneous value of 0min provided that A/V^"(0mm) 
is greater than the Hubble expansion rate. In fact this is the case so long as 
-f^faise > -f^MSSM ~ 1 GcV, scc (|10.22|) . Oucc i^faise ^ -f^MSSM, the flat direction 
potential dominates the energy density of the universe, and MSSM inflation begins 
at a Hubble expansion rate -f^MSSM- In the meantime, landscape tunnelling to vacua 
with smaller cosmological constant continues, and the location of the false minimum 



continuously changes. Eventually \^"(0min) < -^^mssm when 



' false 



2[n — 2) 
at which time 

, , 00 / -f^MSSM 



10.23) 



10.24) 



"""" V. - _ 2)2 V 
It turns out from (|10.17|JT0.18|) that 0min — 0o ^ ^0- Therefore is already inside the 



interval required for a successful MSSM inflation. At this point the Hubble induced 
corrections become largely unimportant, and all the successes of MSSM inflation are 
retained, as discussed in the previous subsection. The fate of the string moduli inside 
the MSSM bubble remains the same as in the previous subsection. In particular, it 



^^Here we are assuming that the difference between the energy densities of the two false vacua is 
small compared to their average energy density. 

■^^Note that tunnellings do not affect the Hubble expansion rate anymore since the flat direction 
dominates the energy density of the universe now. 
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does not matter whether the universe tunnels right away to the currently observed 
value of A or not. Inflation dilutes any excitations of the modulus oscillations during 
inflation. Our Hubble patch reheats when the MSSM fiat direction rolls down to its 
minimum and starts creating MSSM quanta as discussed in the previous subsection. 

Note that all needed for the success of this scenario is to start in a false vacuum 
in the landscape where Hf^ise > rn^ ~ 1 TeV. This ensures that the fiat direction 
will roll way from the origin and settle at 0min, which is the true minimum of V{(j)) 
at that time. It will then track; (pmin as iffaise slowly decreases, and will eventually 
land inside the appropriate interval around (pQ. This is a much milder condition than 
that in the case of negligible Hubble induced corrections -f^faise > 10^ GeV (see the 
previous subsection). 

One comment is in order. This scenario has some similarities to the Affleck-Dine 
baryogenesis with negative Hubble induced corrections [B3, E^. It can be seen from 



( |10.8|) that the equation of motion of the flat direction has a flxed point so long as 
Hubble induced corrections are dominant. The flat direction tracks the flxed point 
in a radiation or matter dominated universe. However, a non-adiabatic change in 
the potential occurs when H ~ m^, and the flat direction starts oscillating around 
the origin. However, in our case the universe is in a de Sitter phase with a slowly 
varying if false, and the flat direction tracks the false minimum of its potential until 
it dominates the universe. As explained, this is necessary for having a successful 
MSSM inflation. 
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